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The thermal decomposition of gaseous aluminum tri- 
methyl, in the absence of hydrogen, is neafly 94 percent 
homogeneous and of kinetic order three-halves over the 
pressure range 10 to 85 mm measured at room tempera- 
ture. The reaction products consist largely of methane and 
to a lesser extent ethane, ethylene, and hydrogen and a 
solid deposit on the walls of the reaction vessel. A reaction 
mechanism involving methyl radicals is proposed. The 
apparent energy of activation, calculated on this basis, is 
found to be 45 kcal., in agreement with the experimentally 
determined figure. The high ratio of methane produced to 
alkyl decomposed suggests that aluminum trimethyl in 


the vapor state at room temperatures must consist at 
least of tetrameric molecular complexes. In the presence 
of hydrogen, both the decomposition rate-and the energy 
of activation are reduced as compared to the decom- 
position of the pure alkyl, a puzzling effect for which no 
satisfactory explanation can at present be offered. Photo- 
chemically, the presence of hydrogen seems to have little 
effect on the over-all decomposition but enhances the 
relative amount of methane slightly. In this, aluminum 
trimethyl appears to differ from mercury dimethyl which 
is known to decompose faster thermally and photo- 
chemically in the presence of hydrogen than in its absence. 





INCE the recent! publication of more precise 

data on the physical properties of aluminum 
trimethyl, such as vapor pressure, molecular 
weight, and molecular association, attention of 
research workers has been drawn to the study of 
the compound, particularly with regard to the 
structure of its fairly stable dimeric molecules.? 
No kinetic study, however, of gaseous aluminum 
trimethyl, has, up to date, been published. It was 
thought to be of interest to investigate the 


* Abstract from the thesis presented by C. C. Schubert 
in partial fulfillment of the requirements for the Master’s 
degree. 
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kinetics of its decomposition and compare the 
behavior of this trivalent metal alkyl with that 
of divalent mercury dimethyl* and tetravalent 
alkyls like silicon tetramethyl‘ and tin tetra- 
methyl,> and obtain at the same time some 
information which would throw light on the 
molecular structure of polymeric aluminum tri- 
methyl, which is still in question. Investigation 
has, therefore, been made of the thermal and 
photochemical decomposition of the compound 
in the gas phase in absence and in presence of 
hydrogen. The kinetic behavior of aluminum 
alkyl appears to be, in many respects, dissimilar 
to that of the corresponding divalent and tetra- 
valent alkyls. 


3 John P. Cunningham and H. S. Taylor, J. Chem. Phys. 
6, 357 (1938). 

4 Helm and Mack, J. Am. Chem. Soc. 59, 60 (1937). 

5 Chas. E. Waring and William S. Horton, J. Am. Chem. 
Soc. 67, 540 (1945). 
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EXPERIMENTAL 
Materials 


Aluminum trimethyl] was prepared® by heating 
together aluminum metal foil and mercury 
dimethyl in a sealed tube on a boiling water 
bath, the mercury alkyl having been obtained 
according to the method indicated by Marvel 
and Gould.’ Since the aluminum alkyl is highly 
and spontaneously inflammable in air, decom- 
posing even in traces of oxygen, it was distilled in 
an oxygen-free nitrogen atmosphere in an all-glass 
set-up. The aluminum alkyl fraction collected 
boiled at 126-127°C at atmospheric pressure. It 
was further purified by distillation in vacuum 
in situ in the gas kinetic apparatus and stored 
in a bulb coated with carbon black and kept in 
a bath of acetone and dry-ice mixture.® 

Hydrogen from commercial cylinders was 
purified by passage over platinized asbestos, 
then through a phosphorus pentoxide tube, and 
finally through a trap cooled in liquid air. A 
quartz mercy lamp of the Uviarc type, run hot 
so as completely to reverse the mercury resonance 
lines, served as the source of illumination in 
photochemical experiments. An electrical furnace 


6 Buckton and Odling, Ann. Chim. Phys. 4, 492 (1865). 

7 Marvel and Gould, J. Am. Chem. Soc. 44, 153 (1922). 

8 Mention may be made of a new alkyl of aluminum and 
mercury which was by accident prepared in this con- 
nection. When the sealed tube containing aluminum foil 
and mercury dimethyl was heated in the water bath for 
several hours beyond the stage at which most of the 
mercury in the alkyl was displaced by aluminum, most of 
the resulting liquid, on fractionation, distilled off * at 
49-51°C at atmospheric pressure, This liquid, from pre- 
liminary tests, appears to be, not a molecular mixture, but 
a true molecular compound of aluminum-mercury alkyl 
and is more highly and spontaneously inflammable than 
the pure aluminum trimethyl. Further study of it is in 


progress. 
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of block aluminum core provided with a quartz 
window for illumination was used for maintaining 
temperature. In the thermal runs, various con- 
stant boiling baths of organic compounds were 
employed. 


Apparatus and Procedure 


The essential features of the apparatus for the 
handling of aluminum alkyl vapor are indicated 
in Fig. 1. The F’s are ground-glass float valves 
by means of which not only could the compound 
be stored in the bulb A out of contact with 
mercury of manometer M, but also prevented 
from contamination by the grease of the stop- 
cocks of the vacuum system. 

The reaction vessel for the thermal runs con- 
sisted of a cylindrical Pyrex tube of 1.9 cm inner 
diameter and 22 cm in length, with a surface to 
volume ratio of about 2. When initially sealed 
at one end’to the manometric system, the reac- 
tion vessel carried at its free extremity a capillary 
inner seal, as shown in Fig. 1. The volume of the 
reaction vessel, and that of any portion of the 
dead space of the manometric system whose 
value was required to be known in order to cal- 
culate the absolute amount of the alkyl in the 
reaction chamber, could be evaluated by ad- 
mitting air or nitrogen into the system and 
measuring the pressures with the help of the 
calibrated McLeod gauge. 

To fill the reaction tube with the desired 
amount of the aluminum alkyl, the storage bulb 
A was maintained at a temperature below room 
temperature ‘so as to avoid coridensation of the 
vapors during the subsequent pressure measure- 
ments, and the vapor admitted into the space 
confined between F), Fs, P, and the reaction 
vessel R, and its pressure was read on the ma- 
nommeter M3. It was then condensed completely 
into R with the help of an ether-dry-ice mixture, 
and the reaction vessel was sealed off from the 
system at S. Knowing the volume of the alkyl 
vapor whose pressure was read on the manom- 
eter and that of the reaction vessel, the pres- 
sure of the vapor in the vessel could be cal- 
culated. When hydrogen was to be added initially 
to the alkyl, the reaction tube was filled first 
with the latter and then with the required 
amount of hydrogen. Similar procedure was 
adopted in the case of the quartz reaction tube 
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which was provided at its two ends with graded 
seals. After the completion of the decomposition, 
the reaction tube was re-sealed to the mano- 
metric system, at the end carrying the inner 
capillary seal which could be broken by an 
iron-in-glass plunger propelled by an electro- 
magnet, in order to study the nature of the 
reaction products. 

It may be noted that before breaking the 
capillary seal as well as before initially filling the 
reaction vessel with the alkyl, the whole mano- 
metric system was flushed with some aluminum 
trimethyl vapor and then completely evacuated. 
This precaution was found necessary in the 
light of preliminary experiments which revealed 
decomposition, to varying degree, of the alkyl 
vapor at first contact with clean glass walls. The 
reaction tubes before use were cleaned, always in 
a uniform manner, with warm chromic acid 
solution, followed by distilled water. 


Analysis 


The reaction vessel containing the products 
of the thermal or photochemical decomposition 
was first surrounded with a liquid-air bath; the 
mixture of non-condensable methane and hy- 
drogen was transferred by means of the Toepler 
T into the McLeod gaugé and its absolute total 
amount measured. The liquid-air bath was then 
replaced by an ether-dry-ice mixture in order to 
remove the two-carbon hydrocarbons and meas- 
ure their amount after transference into the 
McLeod gauge. Since the amounts of the gases 
were extremely small, particularly in the case of 
the two-carbon hydrocarbons at low percentage 
decompositions, their analyses were conducted in 
the microchemical gas analysis apparatus of 
Blacet and Leighton.® In the hydrogen-methane 
mixture, hydrogen was determined by oxidation 
over an alkaline copper oxide bead and the 
methane calculated by difference. In the hydro- 
carbon mixture, ethylene was estimated by 
absorption by fuming sulfuric acid in a porous 
glass bead and the ethane evaluated by differ- 
ence. Tests in this mixture for acetylene by means 
of an alkaline cuprous chloride bead did not give 
positive results. 

°F. E. Blacet, P. A. Leighton, and G. D. MacDonald, 


Ind. Eng. Chem. Anal. Ed. 3, 266 (1931); 5, 272 (1933); 
6, 334 (1934). 


Experimental Results 


The analytical data on the thermal and photo- 
chemical decomposition of aluminum trimethyl 
both in absence and in presence of hydrogen are 
presented in Table I. The pressures of gases are 
given in mm Hg at room temperature. The reac- 
tion time in all cases was 60 min., except in a 
few instances as indicated in the footnotes to the 
table. The last column gives the amount of 
methane produced per unit of alkyl decomposed. 
The zero values represent in some cases traces of 
the particular component. The bracketed experi- 
ments are meant to show the reproducible 
character of a run, which from other data, not 
given in the present table, was quite satisfactory. 
Average values of these experiments will later be 
utilized for calculating specific reaction con- 
stants. Experiments 40 and 39 give the data for 
the decomposition in packed reaction vessels 
with increased surface to volume ratios. The 
unpacked vessels used in all the other experi- 
ments had a ratio of about 2. 

In every run, besides the gaseous reaction 
products recorded in the table, there was a solid 
deposit on the walls of the reaction vessel which 
was in some cases a thin film with a bluish 
metallic sheen, soluble in dilute acid and alkali, 











TABLE I. 

Init. Init. Final CHs 

alkyl He alkyl CHa CoHe CoH, Hz: -—— 
Expt. T7°C mm mm mm mm mm mm mm Alkyl 

Thermal 
40 298 14.7 12.7 9.16 0.69 0.06 0.34 4.6 
39 298 14.9 12.6 12.55 0.92 0.10 0.40 5.4 
38 298 9.1 8.4 3.60 049 001 — 5.1 
32 298 15.4 13.5 6.72 0.57 0.02 0.19 3.6 
34 298 14.8 13.1 685 051 —— 0.21 40 
66¢ 298 15.7 14.4 4.93 049 0.02 0.03 3.8 
6724 298 15.2 0.0 61.70 0.55 0.00 0.00 4.1 
37 298 29.2 24.6 1660 — -— 0.20 3.6 
45 298 29.8 25.0 17.36 0.39 0.02 0.54 3.6 
47 298 48.0 37.8 40.30 0.27 000 —— 40 
48 298 50.9 39.8 43.50 0.43 0.00 0.52 3.9 
49 298 85.5 63.1 205.20 0.60 0.00 0.00 9.2 
, oe | 10.9 17.40 0.52 0.02 0.99 4.3 
65¢ 334 15.1 9.1 27.20 0.58 002 —— 4.5 
$1 334 15.1 4.8 46.20 0.44 0.01 1.10 4.5 
57 298 15.2 148 13.8 8.13 0.32 001 15.2 5.8 
62 298 «615.3 «6155 «64139 8.90 0.33 O01 15.0 74 
71 313 (14.7 15.5 12.6 10.70 O.51 OO1 16.4 5.1 
58 334 14.7 15.1 9.8 28.00 0.38 O01 13.9 $.7 
59 334 148 26.6 9.6 29.40 0.40 0.02 26.4 5.7 
Photochemical 

68¢ 29 «13.8 13.1 0.52 1.47 0.00 0.03 
61 120 13.6 13.1 0.76 0.98 0.00 0.03 
60/ 120 Be 685.7 = =«6(18.8 1.26 1.06 0.00 15.00 








@ In a packed vessel with surface to volume ratio 5.3; > 10.8. 

¢ Decomposition time—40 minutes. 

4 Complete decomposition—several hours. 

¢ Decomposition time—180 minutes; /90 minutes, current being 4 
amp. 
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and in cases of higher temperature and pressure 
and longer reaction period, a complex of metallic 
and brownish polymeric material. At the highest 
pressure investigated, 85.5 mm, a heavy yellow- 
ish deposit was formed at the two ends of the 
reaction vessel; it reacted vigorously with water 
yielding gases which seemed to contain, as 
judged by the smell, some ethylene; the main 
component of these gases might be methane, 
since aluminum carbide, yellowish in color, is 
known to decompose with water into methane 
and aluminum oxide. 

It may be remarked that the chief source of 
error in the data was the measurement of the 
alkyl pressure manometrically, which could not 
be accurate to more than 0.5 mm, so that in 
cases of low percentage decompositions, as in 
Experiment 38, a large error could creep into the 
value of alkyl decomposed, obtained as the dif- 
ference between initial and final pressures re- 
corded in the table. It may be well also to add 
that the tediously long procedure involved in the 
completion of a run is partly responsible for the 
lack of more abundant data at different initial 
pressures of alkyl at 334°C and in the presence 
of hydrogen. 


DISCUSSION OF RESULTS 


Thermal Decomposition in Absence 
of Hydrogen 


Surface Effect 


From the amounts of alkyl decomposed, Ex- 
periments 39 and 40, in packed vessels with a 
surface to volume ratio of 5.3 and 10.8, re- 
spectively, calculation, according to the method 
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indicated by Geddes and Mack,!® leads to about 
6 percent heterogeneity for the reactions con- 
ducted in unpacked vessels employed in our 
experiments whose surface to volume ratio was 
about 2. Consequently the thermal decom- 
position of aluminum alkyl is mainly a homoge- 
neous reaction. 


Order of Reaction 


Preliminary calculations for the order of the 
alkyl decomposition at varying initial pressures, 
employing the expression 


Ax;/Ax2=(p1/p2)", 


where Ax, and Ax, are the decrease in pressure 
per unit time and ; and p»2 are the average 
of initial and final pressures in the two cases, 
gave an average value of about 1.5 for m, the 
order of the reaction. On the basis of three- 
halves kinetic order, the specific reaction rate 
constants for various initial pressures at 298°C 
were calculated using the rate expression in its 
integrated form, 


2 1 1 
k= | ——-<| 
tL(a—x)! a} 


where a and (a—x) are the initial and final con- 
centration, after time ¢ sec., concentration being 
expressed in moles per cc, so that k has the 
dimensions of cc! mole} sec.~!. The constants 
are presented in Table II at three temperatures. 
The slope of the plot of log k versus reciprocal of 
absolute temperature 7, in Fig. 2, yields, for the 
apparent activation energy E of the process, a 
value of about 45 kcal. The frequency factor A, 
in the rate expression, 


k=A Xe~15000/RT 


has a value of 8.1X10'*, quite comparable as to 





TABLE II. Percentage composition of gaseous products at 
334°C with an initial alkyl pressure of 15 mm. 











Time in 

Expt. min. CH, CoHe C2Ha He 
53 5 91.34 8.03 0.63 Traces 
54 10 93.36 6.16 0.48 Traces 
52 15 95.07 4.80 0.12 Traces 
65 40 97.84 2.09 0.07 Traces 

51 60 97.00 0.90 0.04 1.9 

9 180 98.00 0.50 0.00 9 








10 Geddes and Mack, J. Am. Chem. Soc. 52, 4372 (1930). 
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the order of magnitude with the values of fre- 
quency factors for other three-halves order reac- 
tions." The slight increase over the normal value 
of 10'5 is not unexpected in view of the probably 
high entropy of activation which can be ac- 
counted for by the polymeric nature of the 
reactant in its initial state, as will be presently 
indicated. The three-halves order might, at 
first sight, appear to hold only within the pres- 
sure range recorded in Table I, namely 10 to 
85.5 mm. Actually, however, since aluminum 
trimethyl is known to exist completely in the 
dimeric form even at 70°C,! and since the experi- 
mental pressures were measured at room tem- 


perature at which the alkyl might be in a higher - 


polymeric state, the effective pressure range, in 
reference to the decomposition temperatures at 
which the polymer molecules would be com- 
pletely dissociated into simple monomeric mole- 
cules, would be about 20 to 170 mm. Over this 
whole range, therefore, the kinetic order of the 
reaction would seem to be three-halves. 


Reaction Products 


The predominant component of gaseous reac- 
tion products is methane, with smaller amounts 
of hydrogen, ethane, and ethylene, besides a 
complex solid deposit on the walls of reaction 
vessel, which might be composed of metallic 
aluminum, aluminum carbide, and some polymer 
material. The variation of the percentage com- 
position of gaseous products with time at 334°C 
is expressed by the data in Table III. Impossi- 
bility of obtaining measurable and analyzable 
amounts of products prevented the performance 
of runs with shorter periods than 5 min., and at 
temperatures below 334°C. The general trend 
of the results obtained is nevertheless quite clear. 
Methane percentage, and in a smaller degree 
that of hydrogen, increase with time, whereas the 
percentages of ethane and ethylene suffer a 
regular decrease. This may be attributed to 
progressive pyrolysis with time of ethane to give 
methane, thus increasing total methane per- 
centage. The gradual decrease of ethylene and 
its final disappearance in three hours, Experi- 
ment 9, might be due, partly, to its hydrogena- 


11 See R. N. Pease, Equilibrium and Kinetics of Gas Reac- 
tions (Princeton University Press, Princeton, 1942), p. 133. 


TABLE III. 








T°C 298 298 298 298 313 334 


Init. ' 
alkyl 15.1 29.5 49.5 85.5 15.1 15.1 


mm 
kis 0.042 0.040 0.044 0.042 0.110 0.476 








tion, and partly to pyrolysis and subsequent 


‘ polymerization. 


It is of interest to note the high ratio of 
methane to alkyl, about 9, in the last column of 
Table I, at the highest initial pressure of alkyl 
employed, namely 85.5 mm. The repetition of 
another high initial pressure run gave also a 
similar high methane to alkyl ratio. This pre- 
supposes at least a tetrameric form of the 
aluminum alkyl in the vapor at room tempera- 
ture, as will be evident from the following 
stoichiometric relation, 


[Al (CH 3) 3 Ja —AlL,Cs +9CH 4. 


The possibility of such a relatively high poly- 
meric form of aluminum trimethyl is suggested 
by the nature of the molecular weight versus 
temperature curve of Laubengayer and Gil- 
liam,! and also by the ability of this alkyl to 
form associated complexes with other com- 
pounds.” Hence aluminum trimethyl seems to 
exist in the vapor phase at room temperatures in 
a more complex molecular state than dimeric, 
probably in a tetrameric form. 


Reaction Mechanism 


To account for the above facts, the following 
mechanism is suggested. The thermal decom- 
position of aluminum trimethyl being of three- 
halves order, must, by analogy with other similar 
reactions, involve a chain mechanism, the initial 
step being the splitting of a methyl radical 
according to the equation, 


Al(CH3)3—Al(CH3)2+CHs. (1) 


The complex radical Al(CHs3)2 might break down 
further into methyl radicals or alternately could 
be the source of ethane which appears in rela- 
tively high percentages in the early stages of 
the reaction, before it has time to decompose as 
in the long period runs in which it is present in 
smaller percentages. This would also account for 


12 Norman Davidson and Herbert C. Brown, J. Am. 
Chem. Soc. 64, 316 (1942). 
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the metallic deposit observed on the walls of the 
reaction vessel. 

The formation of methane probably occurs by 
the reaction, 


Al(CH3)3+CH;—~CH,+Al(CHs3)2:CHs. (2a) 


The latter radical, by analogy with other com- 
plex radicals, like propyl, etc., may tend to lose 
a methyl radical and revert to a stable molecule. 


Al(CH;)2:CH:>CH;+Al(CH;)-CH». (2b) 


Actually compounds of the type BrAl: CH2 have 
been prepared from aluminum and CH2Bryz. It 
seems probable therefore that Al(CH;)-CHe 
could be stable and could deposit out in this or 
in a polymerized form on the walls of the reaction 
vessel, thus accounting for some of the solid 
deposit which has been remarked earlier. 

Steps (2a) and (2b) would together constitute 
a chain. The chain breaking step must be the 
recombination of the radicals into ethane, 


CH;+CH;—C2H,g, (3) 


in order to account for three-halves order. 
Further reactions of the type (2a) and (2b). with 
the eventual production of aluminum carbide, 
the presence of which in the deposit was indicated 
by its vigorous reaction with water, as has 
already been remarked, would kinetically be 
similar. The molecule in step (20) in its simple 
or polymeric form would probably also react 
with water to yield some methane and plausibly 
also unsaturated hydrocarbons. For all these 
reactions a constant k, may be assigned, which 
would evidently. be a complex of the rates of 
several reactions of the type (2a). 

The steady-state concentration of methyl 
radicals, on this basis, is given by steps (1) and 
(3), as the steps (2a) and (20) mutually cancel. 


[cH,]= (191) 


[A] standing for aluminum alkyl concentration. 
The over-all. rate of disappearance of the alkyl 
is given by 
a[A] 
——a = LA] +h[A (CH). 
Substitution for concentration of methyl radicals 


yields 
a[A] 


~ Sh -nta]tn(2) cA? 
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The initial rate k,[A] will be negligible in com- 
parison with the rate of chain reactions, and 
therefore the observed three-halves order for the 
over-all decomposition is accounted for. 

On this basis, the apparent energy of activa- 
tion E observed will be 


E=E.+3(Ei—£5). 


E, the energy of activation of reaction (2a) is 
that of a radical-molecule reaction, which seems 
to have in general a value of 15 kcal. or less. 
E;, the energy of activation of recombination of 
two methyl radicals is very small and may be 
neglected. E,, essentially the strength of the 
Al—C bond in the alkyl is not known. A probable 
estimation can be made from a knowledge of 
related alkyl compounds. For the Hg—C bond 
in mercury dimethyl, Terenin™ gives a value of 
40-50 kcal. Sodium methyl decomposes rapidly 
at 100°C, and mercury dimethyl would decom- 
pose rapidly around 400°C.’ This would mean a 
value of 20 to 25 kcal. for the Na—C bond. The 
C—C bond energy is about 85 kcal. A value of 
about 60 kcal. for the bond between C and an 
element like Al in group III of the Periodic 
Table seems reasonable. Substitution of this 
value for E; in the expression will give for E of 
the over-all decomposition 


E=15+43(60)=45 kcal. 


a value in agreement with the experimental value 
of ~45 kcal. 

Compared to tin and silicon tetramethyls, 
aluminum ‘trimethyl appears to be different as 
regards both kinetic order and apparent energy 
of activation. Waring and Horton® report for 
thermal decomposition of tin tetramethyl a 
first-order constant above 80 mm of initial alkyl 
pressure and below this pressure an order ap- 
proaching 2, and an energy of activation for the 
first-order reaction of 82 kcal., which high value 
is interpreted as evidence of a decomposition 
mechanism through molecular rearrangement. 
Silicon tetramethyl‘ too decomposes according 
to first order, with an energy of activation of 
79 kcal. For dimethyl mercury,’ similar data 
are not available; but it is known that it decom- 
poses readily only at a higher temperature than 


13 A. Terenin and N. Prilezhaeva, Acta Physicochimica 
1, 759 (1934). 
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does aluminum trimethyl. An adequate explana- 
tion of these differences would require further 
research and examination 


Thermal Decomposition in Presence of 
Hydrogen 


The effect of hydrogen, as is seen from the 
last section of the data in Table I, is to decrease 
the over-all rate of decomposition of the alkyl 
and to increase the ratio of methane to alkyl 
decomposed. An approximate doubling of hy- 
drogen pressure for the same initial alkyl pres- 
sure of about 14.7 mm, cf. Experiment 59, 
produces little change over the value of Ex- 
periment 58. With this piece of evidence, in 
spite of the absence of data on the effect of 
varying initial pressure of alkyl in presence of 
a constant pressure of hydrogen, it might ten- 
tatively be assumed that the kinetic order of 
alkyl decomposition is three-halves as in the 
case of the pure alkyl. The specific reaction rate 
constants calculated on this basis are presented 
in Table IV at three temperatures, and the plot 
of log k versus reciprocal of absolute temperature 
is given in Fig. 3. The slope of this straight line 
yields for the apparent energy of activation E 
a value of 30.5 kcal. Calculation of E on the 
basis of a first-order reaction, just as a check, 
gave a value of 29 kcal., not materialiy different 
from the previous one. This is rather puzzling in 
that the presence of hydrogen should reduce 
both the rate of decomposition and the apparent 
energy of activation. At present we are not able 
to offer any satisfactory explanation. The fact 
that the presence of hydrogen increases the ratio 
of methane to alkyl decomposed from about 4 
to 5, does not throw much light on the above 
mentioned problem. It might well be that the 
kinetics in the two cases are quite different and 
can be elucidated only by further investigation 
of the problem. 

It must be noted that the behavior of alu- 
minum alkyl in presence of hydrogen is markedly 








TABLE IV. 
: 298 313 334 
Init. alkyl mm 15.1 14.7 14.7 
Ris 0.028 0.050 0.140 











+S) 


Log k . 








1 i ! 
162 166 170 Lm 
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different from that of mercury dimethyl* which 
decomposes faster in presence of hydrogen than 
in its absence, the accelerated rate increasing 
with increase in hydrogen pressure and con- 
suming hydrogen in the process. In the case of 
aluminum, the adverse effect of hydrogen be- 
comes constant after a pressure value of 15 mm, 
further increase of hydrogen making little dif- 
ference, as can be seen from Experiments 58 and 
59, though no hydrogen is used up. 


Photochemical Decomposition 


The few photochemical experiments whose 
results are presented in Table I were performed 
primarily to verify whether the difference which 
was observed in the thermal decomposition of 
aluminum and mercury alkyls would hold also 
in the photodecomposition. The effect of hy- 
drogen in Experiment 60 as compared with its 
absence in Experiment 61, is not appreciable at 
all, whereas in the case of dimethyl mercury, 
hydrogen noticeably accelerated the rate of de- 
composition of the alkyl. No definite conclusions 
can now be drawn regarding the nature of the 
photochemical decomposition of aluminum alkyl 
in absence of exhaustive data. Still the difference 
in behavior of aluminum and mercury alkyls is 
evident. Further photochemical study of alu- 
minum trimethyl will soon be undertaken in 
this laboratory. 
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A method is described that permits deduction of selection rules for Raman and infra-red 
spectra of crystals from a knowledge of just the space group designation and the population of 
the unit cell. When implemented by an appendix to this article showing, in a convenient 
arrangement, the distribution of point symmetries within space groups the method can be 
applied by a user familiar with point group manipulation but having no special knowledge of 
space groups. There are no selection rules operating in the liquid state. Some tentative qualita- 
tive conclusions are drawn concerning the relative intensities of certain kinds of components 
appearing in spectra of condensed systems. Sharpened selection rules are expected, and have 
been observed by others, to operate for single crystals containing non-rotating molecules when 
these are given special orientations in the light path. These special orientations and the appro- 
priate selection rules can be predicted from information mentioned above along with inspection 
of the external form of a crystal. All procedures have a desirable advantage of simplicity over 
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prior ones directed toward the same ends. 





INTRODUCTION 


ETAILED studies, both experimental and 
theoretical, of the motions of isolated mole- 

cules have yielded a vast dividend of valuable 
accurate information about gaseous substances. 
This information not only has immediate physical 
significance but also, when combined with power- 
ful tools of statistical mechanics, contains many 
chemical implications. The store of substances 
which may be regarded conveniently as vapors is 
becoming rapidly exhausted while generally re- 
liable procedures for dealing similarly with sub- 
stances in condensed states remain to be de- 
veloped. Theoretical problems: which must be 
overcome in this connection are severely compli- 
cated and are not likely to be solved without the 
aid of simplifying idealizations. The possibilities 
in this direction are numerous while an intelligent 
assessment of their merits is hampered by an 
insufficiency of experimental information in 
several categories, especially concerning the de- 
tails of molecular motions in condensed systems. 
The most direct experimental route to knowl- 
edge of molecular motions in any situation is 
through the spectra they produce. The interpre- 
tation of spectra demands certain concomitant 
theoretical machinery. This communication will 
be devoted to the development of some principles 
appropriate for condensed phases, that will be 
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applied subsequently in separate reports to the 
interpretation of experimental investigations. 


SELECTION RULES FOR CRYSTALS 


Well-known procedures! permit the deduction 
of selection rules for vibration-rotation spectra of 
isolated symmetric molecules. To review them 
briefly, the mutual potential energy of a collec- 
tion of interacting atoms will reflect the sym- 
metry of their equilibrium configuration. This 
symmetry can be described by means of a group 
of physical operations. The term group is used 
here in a specialized mathematical sense and 
implies many consequences that follow logically 
from the restrictive definition. For instance, 
associated with the group is a set of irreducible 
representations which may be visualized as non- 
equivalent examples of the given symmetry. 
These are limited in number and have the 
property that all other examples of the same sym- 
metry can be constructed from combinations of 
the irreducible ones. Thus, the potential energy 
of a distribution of interacting points can be 
analyzed into a sum of independent contributions. 
Each contribution is associated with some one of 
the irreducible representations of the group 
describing the symmetry of the equilibrium con- 


1 For details see J. E. Rosenthal and G. M. Murphy, Rev. 


=e (1936); L. Tisza, Zeits. f. Physik 82, 48 
1933). 





mpre errerces 


ee ee ee eee 


eT re as 


Peas 


See Sas 5 FO 





SLAG ONS RS ES 








S, 





Be BL, Sia ace: set 


EAT BETES 





INTENSITIES IN CONDENSED PHASES 9 


figuration. Each contribution to the analyzed 
potential energy determines a set of fundamental 
modes of motion. The totality of fundamental 
modes is the composite sum extended over all 
irreducible representations. Finally, the collection 
of points in motion will possess a transitory elec- 
tric moment and a transitory polarizability. 
These properties furnish two mechanisms for 


exchange of energy with a radiation field and so , 


give rise to infra-red and Raman spectra, re- 
spectively. These properties also reflect the sym- 
metry of the equilibrium configuration and there- 
fore can be analyzed into components associated 
with the irreducible representations. Only those. 
modes of motion will be concerned in the corre- 
sponding spectra that appear in irreducible 
representations along with components of the 
transitory electric moment or transitory polariza- 
bility. Selection rules deduced in this way from 
the symmetry of an isolated molecular configura- 
tion find application in the analysis of spegtra of 
dilute vapors where the condition of isolation is 
achieved. However, they are inappropriate for 
the analysis of spectra of condensed phases. 

It is evident from the foregoing discussion that 
selection rules can be deduced in the same way 
for any situation where it is possible to specify 
and analyze the symmetry ofa suitable configura- 
tion of interacting points. Recognizing this op- 
portunity, Bhagavantam and Venkatarayudu’ 
have taken the lead in developing methods to 
permit deduction of selection rules for molecular 
crystals. They propose that all of the spectro- 
scopically important frequencies can be discovered 
and classified by examination of the isolated 
crystallographic unit cell. This proposal is almost 
certainly correct, not only for the reasons ad- 
vanced by its authors but also in the light of 
experimental findings. The unit cell is then 
regarded as an example of the space group ap- 
propriate for the crystal lattice ; a point is deemed 
invariant under any operation that carries it 
either into itself or into the identical point in 
another unit cell but not so under an operation 
that carries it into an equivalent point in its own 
or another unit cell. The logic of this last con- 
vention, considered in relation to the potential 


2S. Bhagavantam and T. Venkatarayudu, Proc. Ind. 
Acad. Sci. 9A, 224 (1939); S. Bhagavantam, Proc. Ind. 
Acad. Sci. 13A, 543 (1941). 


energy of the complete crystal, is not altogether 
satisfactory under close scrutiny. 

This theoretical machinery for crystals has 
been tested against a number of experimental 
cases.2* In practice it is found that the method 
introduces and emphasizes distinctions that are 
without meaning in relation to available experi- 
mental techniques. The results require further 
assessment and interpretation on this account so 
that the method is often cumbersome and some- 
times even tedious to apply. The rhombic modifi- 
cation of crystalline sulfur furnishes an extreme 
example. The unit cell contains 16 molecules of 
Ss and exhibits 384 distinct modes of motion 
whereas the Raman spectrum consists of only 10 
lines. This discrepancy is not the result of ex- 
ceedingly strict selection rules but instead is 
caused by the circumstance that the 384 modes 
can be collected into sets such that the theo- 
retically distinguishable frequencies in any one 
set are indistinguishable in practice. In a mo- 
lecular crystal it is convenient to consider 
separately two types of motions: (1) lattice 
modes, which are approximately the motions of 
rigid molecules and (2) molecular modes, which 
are approximately the distortions of molecules 
whose centers of mass and principal axes of 
inertia are at rest. For sulfur, 288 molecular 
modes are distributed among 14 sets, each set 
having a representative frequency not sensibly 
different from one of the 14 fundamental fre- 
quencies of an isolated Ss molecule. The members 
of a given set correspond to different ways of 
coupling the motions of 16 molecules when all are 
moving with the same frequency. In principle the 
coupling will yield 16 new frequencies but in 
practice the distinctions are not apparent. Similar 
considerations apply to lattice modes although in 
this category the opportunities should be better 
for practical resolution of the distinctions. Even 
so, in the example under discussion, only 2 some- 
what broadened lines can be attributed to 48 
modes arising out of restricted rotation of mole- 
cules contained in the unit cell. The lattice modes 
of translatory origin have low frequencies along 
with low intensities and so are not likely to be 
observed. In general the organization of modes 


( AD) review by G. Glockler: Rev. Mod. Phys. 15, 111 
1943). 
4R. Norris, Proc. Ind. Acad. Sci. 13A, 291 (1941). 
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into sets can be accomplished more or less readily ; 
when this has been done it still remains to ex- 
amine the different selection rules for the mem- 
bers of a set in order to arrive at a representative 
selection rule for the common frequency. In 
summary it may be said that this ingenious 
procedure leaves something to be desired, both 
logically and practically speaking. 

- This writer prefers to employ another plan of 
attack. The idealization adopted is the one so 
strongly implied in experimental findings; i.e., to 
treat the motions of one molecule moving in a 
potential field reflecting the symmetry of the 
surrounding crystal. This permits the use of point 
symmetry rather than space symmetry and 
attributes a more logical quality of invariance to 
the potential energy of the complete crystal 
under the appropriate symmetry operations. The 
treatment yields at once a suitable number of 
frequencies with appropriate selection rules. Of 
course, the corresponding motions are fictitious in 
principle but so are those attributed by the other 
procedure to an isolated unit cell. If one is forced 
to adopt an idealization there is little to be gained 


by making it unduly elaborate. Furthermore, 


there is good reason to believe that the picture of 
molecular modes derived from examination of 
one molecule may be more accurate than the one 
derived from examination of the unit cell; 
coupling between molecular motions in different 
molecules is extremely weak, ex post facto. 

A prerequisite to the treatment of any molecu- 
lar crystal is a suitable knowledge of the crystal 
structure. This is provided by a detailed pattern 
of the unit cell, when one is available, but most 
often when found in convenient summaries con- 
sists merely of the space group designation along 
with the number of molecules per unit cell. It is 
highly expedient to have theoretical machinery 
designed to manufacture the desired results from 
the latter form of raw materials since these are so 
comparatively easy to procure. Such a formalism 
can be developed to operate in the great majority, 
if not all cases. If it is found to fail occasionally, 
the shortcoming will be apparent at once and 
resort can then be had to the information con- 
tained in a detailed pattern of the unit cell. 

In order to make clear the basis for the 
formalism it is convenient to recall at this time 
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a few well-founded principles of crystallography.® 
Crystals can be assigned, on the basis of their 
external symmetries, among 32 classes, each class 
being identified with a collection of symmetry 
elements and a unique point group of operations 
concerned with them. Each point group generates 
a characteristic number of space groups, de- 
scriptive of the internal crystal structure, when 
its operations are combined in various numbers 
and various ways with the primitive translations 
that carry one lattice point into another. An 


_arbitrarily chosen point-in the crystal is a mem- 


ber of a set of equivalent points, all contained 
within the same unit cell. The set can be gener- 
ated by applying in turn the operations of the 
space group to the chosen point. In general, the 
number of equivalent points per set is equal to 
the number of operations in the space group. 
However, if a point is situated on one of the 
symmetry elements and if the operations of the 
space group concerned with that element are not 
composite translations, then the point remains 
invariant under these operations. Consequently, 
such a point belongs to a set having population 
less than the number of operations in the space 
group. Although usually called otherwise, such a 
point will be designated henceforth as a site. In 
general, a unit cell exhibits several different kinds 
of sites and sometimes several distinct sets of the 
same kind. As a further consequence of its 
definition a site is necessarily the affix for a point 
group, the site group, which is a sub-group of the 
space group. Thus, an infinite crystal (the limit 
is always appropriate in practice) has the point 
symmetry of a site group about each and every 
site. 

The center of mass of a molecule is invariant 
under the operations of the associated molecular 
group and its equilibrium position in the crystal 
will be called the affix of this group. It is a well- 
substantiated rule’ of crystallography that the 
affixes of symmetrical molecules are usually 
situated on sites. In order that this be so it is 
necessary that the site group be a sub-group of the 
molecular group. Thus, the complete potential 
energy of a crystal when considered from the affix 


5 For detailed development see R. W. G. Wyckoff, The 
Analytical Expression of the Results of the Theory of Space 
—_— (Carnegie Institution of Washington, 1930), second 
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of any of its constituent molecules always has the 
symmetry of the appropriate site group. When 
the site group has been identified the selection 
rules can be deduced in the usual way for its 
irreducible representations. Furthermore, since 
the occupant of the site is a valid example of its 
symmetry the motions of one molecule can be 
assigned to the irreducible representations with- 
out regard for the motions of the rest of the 
crystal. Since this procedure accounts for a 
suitable number of frequencies no further elabora- 
tion is required. Of course, one retains the 
reservation that each motion of the one molecule 
is coupled to some degree with the motions of the 
rest of the crystal that are associated with the 
same irreducible representation of the site group. 
The relation between motions of an isolated 
molecule and those of one in the crystal can be 
established readily by comparing character tables 
for the molecular group and the site group. 

The deduction of selection rules for a molecular 
crystal can be accomplished as outlined above, 
whenever the site groups have been identified. An 
acceptable site group must be a sub-group of both 
the space group and the molecular group. When 
the latter two groups are compared only a few 
common sub-groups will emerge and sometimes 
only one. In the latter case the choice is unequivo- 
cal and the site group is unambiguously identified. 
When more than one acceptable site group 
emerges, a further elimination is possible on the 
basis of the number of molecules contained in the 
unit cell. The molecular affixes must occupy a 
complete set of sites and the number per set 


‘ varies with the order of the site group. When 


applying this criterion it is necessary to take into 
account the possible multiplicity of identical sets 
of sites, for several such sets may be in use. In 
most cases the foregoing considerations will 
suffice, but ambiguous situations are possible. 
These arise infrequently; when they do it is 
perhaps safest to refer to a map of the unit cell. 
In order to employ the easy procedure indicated 
in this paragraph it is necessary to have at hand a 
tabulation of the site groups contained in each 
space group, the numbers of equivalent sites per 
set, and the multiplicities of sets. This informa- 
tion can always be found in works on crystal- 
lography, but not arranged conveniently for the 
present purpose. Therefore, in order to implement 


the procedure to the fullest possible extent, the 
requisite material has been compiled in a con- 
venient form and appears as an appendix to this 
article. 

An analysis of the spectra of calcite and 
aragonite will be undertaken now as an illustra- 
tion of the procedure. Calcite is a representation 
of the space group D;,° with two molecules per 
unit cell. Consultation of the appendix shows 
that this population can be accommodated only 
on site groups D; and C;;. The first is a sub-group 
of the molecular group Dy, for carbonate ion 
whereas the second is not. On the other hand the 
second is suited to the point symmetry of calcium 
ion. Hence carbonate ion is surely on site D; and 
calcium ion on site C3;;. Resorting next to the 
usual analysis of group characters! one finds that 
the three lattice modes of translatory origin 
associated with calcium ion are forbidden in the 
Raman spectrum but permitted in the infra-red 
spectrum. Similarly one finds that carbonate ion 
yields four molecular frequencies, one permitted 
only in the Raman spectrum, another permitted 
only in the infra-red spectrum, and two others 
permitted in both spectra. In addition it yields 
one lattice frequency of translatory origin and 
one of rotatory origin active in both spectra along 
with one frequency of each type active only in the 
infra-red spectrum. The conclusions concerning 
molecular frequencies are in complete harmony 
with experimental observations while the ones 
concerning lattice frequencies are not violated 
although fewer than the expected number of 
frequencies are observed.” The contrasting con- 
clusions for aragonite follow readily from the 
knowledge that this crystal is a representation of 
the space group V,!° with four molecules per unit 
cell. Reference to the appendix places calcium ion 
on site C; and carbonate ion on site C,. The 
lattice modes attributable to calcium ion are once 
again forbidden in the Raman spectrum but 
permitted in the infra-red spectrum. However, 
six molecular frequencies are now permitted to 
appear in both spectra. The two extra frequencies 
arise from removal of a degeneracy prevailing in 
both the free carbonate ion and calcite; if the 
splitting is not resolved only four frequencies will 
be observed. The lattice modes attributable to 
carbonate ion are now six in number and are 
permitted to appear in both spectra. Once more, 
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these conclusions are in harmony with experi- 
mental data.? Essentially the same predictions 
can be made, but less conveniently, by the 
method employing space groups. 

The only failure of the easy procedure, en- 
countered thus far by the writer in the course of 
examination of a substantial number of crystals, 
occurs with rhombic sulfur. A representation of 
the space group V,”, the unit cell, contains 
sixteen molecules. Referring to the appendix one 
finds that this population can be accommodated 
equally well on two sets of eight sites each with 
symmetry V or on one set of sites with symmetry 
C;. Further discussion will be omitted except to 
mention that the analysis now involves only 24 
instead of 384 modes. 

In general, the site group will be of lower order 
than the molecular group although this is not 
necessarily the case. When true, the selection 
rules for the crystal will be less strict than those 
for the isolated molecule with the result that the 
crystal will possess a richer spectrum than the 
vapor. For site groups C2, C,, and C3, it can be 
said without further examination that no selec- 
tion rules will operate. 

An interesting aberration develops in the case 
of ammonium chloride. The ammonium ion 
occupies a site having point symmetry O, which 
is of higher order than the molecular group T.. 
From this circumstance one must infer that the 
ion interacts with the rest of the crystal as though 
it were a point, a view which is entirely consistent 
with the prevalent conception that the ion is able 
to rotate freely in this crystal. Special selection 
rules will apply. 


LIQUIDS 


Turning now to the other common condensed 
state, there will be no selection rules operating in 
the liquid phase. This conclusion follows at once 
from the conception of the liquid state as one of 
disorder. The short range order in a liquid does 
not connote symmetry. Convincing proof will be 
presented in a separate communication describing 
an experimental study of liquid and solid benzene. 
For the present, attention is directed toward the 
observation® that the infra-red spectrum of 

* Private communication from Dr. Daniel Norman, New 


England Spectrochemical Laboratories, Medway, Massa- 
chusetts. 
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pentaerythritol, C(CH2OH),, in solution exhibits 
a characteristic absorption band in a certain 
neighborhood whereas a comparable amount of 
solid in the light path is transparent at the same 
place while exhibiting other regions of absorption 
in common with the solution. Molecules of 
pentaerythritol’ occupy sites S, in the crystal. 
This group has four irreducible representations 
and one of them, totally symmetric, containing 
six frequencies is forbidden to contribute to the 
infra-red spectrum. The case does not afford 
unequivocal proof of the general absence of selec- 
tion rules for liquids. It can be argued that 
hydrogen bonds stabilize a symmetrical form in 
the crystal whereas the molecule, isolated from 
others of its own kind, might exhibit several con- 
figurations of approximately equivalent energy 
and no net symmetry. This alternative seems 
unlikely when attention is given to the observa- 
tion® that pentaerythritol halogenohydrins have 
point symmetry S, in the vapor phase. 


RELATIVE INTENSITIES 


Although selection rules forbid certain modes 
of motion from contributing to a given spectrum 
they do not guarantee that contributions from 
the others will always be observed. It is outside 
the scope of currently workable theory to predict 
quantitatively the relative intensities in different 
regions of a spectrum. Nevertheless, cogent argu- 
ments can be directed at some aspects of this 
problem for a substance which follows the domi- 
nant pattern of behavior. Certain modes of 
motion develop spectroscopic activity in a con- 
densed phase but do not contribute to the 
spectrum of the vapor. The intensity of this newly 
acquired activity depends upon the degree of 
interaction between the molecule and its sur- 
roundings. In a similar manner, for those modes 
of motion that are permitted to contribute to 
both spectra, the intensities depend as well upon 
the interactions residing within the molecule. In 
any condensed system where it is legitimate to 
speak of a collection of molecules the latter class 
of interactions must be considerably stronger than 
the former one. Hence, the contribution to-the 
spectrum of a condensed phase, by a mode that 

7F, J. Llewellyn, E. G. Cox, and T. H. Goodwin, J. 


Chem. Soc. 882 (1937). 
8 H. de Laszlo, Comptes rendus 198, 2235 (1934). 
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is active in the spectrum of the isolated molecule, 
will certainly be more intense than the contri- 
bution from a mode that becomes active only 
after condensation. Thus, the customary practice 
seems quite legitimate whereby the most intense 
frequencies are selected from the spectrum of a 
liquid for assignment to fundamentals permitted 
by the selection rules for the isolated molecule. 
At the same time, the interactions causing con- 
densation are substantial, and it is an established 
fact that new contributions of appreciable in- 
tensity do appear in the spectrum of the con- 
densed phase. Furthermore, the interactions 
become more powerful in a crystal than ina liquid 
so that some general increase of intensity might 
be expected to accompany solidification. This 
increase, when observed, will be most noticeable 
for fundamentals that develop activity only after 
entering the condensed phase. Components of the 
spectrum permitted also in the vapor will exhibit 
merely a small increment in a basic level of 
intensity; combinations that develop only after 
entering the condensed phase will be compara- 
tively weak so that the increase of intensity will 
not be striking; the newly developed funda- 
mentals will have medium intensity and their 
growth will seem more spectacular. Lattice modes 
of rotatory origin can develop high intensities,* 4 
whereas those of translatory origin may be strong 
in the infra-red spectra of ionic crystals but are 
probably weak otherwise. 


ORIENTATION EFFECTS 


The selection rules deduced for crystals by 
procedures described here are sometimes sus- 
ceptible to further refinement for special con- 
ditions of observation. In addition to the usual 
alterations of a Raman spectrum resulting from 
the use of polarized light there are entirely new 
possibilities that have no counterpart in spectra 
of fluids. When a fluid is irradiated, its randomly 
oriented molecules present all possible aspects to 
the radiation. Consequently, the electric vector 
of the light can enter all suitable irreducible 
representations of the molecules in the course of a 
sufficiently large number of molecular encounters. 
However, in a large single crystal all molecules 
not rotating freely have similar orientations so 
that a beam of parallel light encounters much the 
same aspect of all molecules in its path. In some 


cases, depending upon the nature of the sie 
groups involved, a direction of approach can be ° 
selected so that no component of the electric 
vector will appear in one of the suitable irre- 
ducible representations of sites that are occupied. 
Then, only two of the components of transitory 
electric moment and only three of the elements of 
transitory polarizability will be energized. For 
some site groups it is possible in addition to select 
a plane of polarization for the incident light so as 
to exclude its components at the same time from 
another irreducible representation. Then, only 
one component of the transitory electric moment 
and one element of the transitory polarizability 
will be energized. The correspondingly sharpened 
selection rules and the directions of approach and 
planes of polarization for which they operate can 
be deduced easily from inspection of the external 
form of the crystal (provided that it is not 
pseudomorphous) and knowledge of the site 
groups that are in use. The irreducible repre- 
sentations of a site group that can be penetrated 
by a vector representing the radiation field will 
be the ones containing components of the transi- 
tory electric moment. The orientations of these 
components relative to the elements of symmetry 
of the site group can be established by inspecting 
the table of group characters. However, the 
elements of symmetry contained in the site group 
must be identical with those contained in the 
crystallographic point group and hence they can be 
located relative to the faces of the crystal by 
inspecting its external form. To predict the effect 
produced in the Raman spectrum is likewise 
easy. Each element of polarizability is associated 
with two, not necessarily distinct components of 
the electric vector of the radiation field. When 
two irreducible representations, not necessarily 
distinct but each containing a component of the 
electric vector, are combined in a direct product! 
the resulting irreducible representation will con- 
tain the associated element of polarizability. If 
either of the components is absent for a particular 
orientation of the crystal then the associated 
element of polarizability will be inoperable. The 
manipulation of group characters required for 
this analysis is exactly the same one used to 
predict selection rules for combinations and 
overtones and would be undertaken in any event. 
The totally symmetric irreducible representation 
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of a site group will often contain the three 
diagonal elements of the polarizability tensor. 
Then this representation will contribute to the 
. spectrum for all orientations of the single crystal. 
This persistence notwithstanding it is evident 
that relative intensities observed in spectra of 
single crystals will be highly sensitive to orienta- 
tion. The validity of this conclusion is well 
established in fact.**%!© The entire subject has 
received elaborate theoretical treatment.’ It is 
interesting to note that a crystal containing 
freely rotating ions should not show the expected 
effects of orientation. This suggests a valuable 
tool for study of some second-order phase 
transitions. 


FREQUENCY ASSIGNMENTS 


Finally it will be remarked that the most 
fruitful approach, when feasible, to the identifi- 
cation of vibration frequencies for a symmetrical 
molecule will start with the spectra of the vapor, 
proceed to the spectra of the solid, and conclude 
with the spectra of the liquid. In favorable cases 
a complete assignment of frequencies can be 
made in this way without having any resort to 
data for isotopic molecules. 


APPENDIX 


Availabilities and Symmetries of Sites in the 230 
Space Groups® 


The standard symbol for the space group will 
be found below in boldface type. Following it are 
a series of symbols describing the various sites. 
The description for a site consists of the standard 
point group symbol preceded by the number of 
distinct sets of this kind and followed in paren- 
theses by the number of equivalent sites per set. 
The numbers are appended only when they differ 
from unity. Thus, the entry C,!: 2C,, C2(2) 
signifies that the space group C,! contains 2 
distinct sets of sites having point symmetry C, 
with one site per set and 1 set of sites having 
point symmetry C, with two equivalent sites per 
set. If the symbol for a space group is absent this 
means that it contains no sites. 


C2:8C,;C,!:2C,; C,?: C,(2); Co!: 4C2; Co®: 2C2(2); 
C.,!: 8Cx,, 4C,(2), 2C,(2); Co,?: 4C;(2), C,(2); 


9S. Bhagavantam, Proc. Ind. Acad. Sci. 11A, 62 (1940). 
10 R, B. Barnes, Phys. Rev. 39, 562 (1932). 
1 J. Cabannes, Comptes rendus 211, 625, 750 (1940). 
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Con*: 4C2n,(2), 2C:(4), 2C2(4), Co(4); Cont: 4C;(2), 
2C2(2); Con’: 4C:(2); Carx®: 4C;(4), C2(4); Car!: 
4Cov, 4C,(2); Co?: 2C,(2); Co?: 4C2(2); Cav!: 
2C2(2), C.(2); CoS: 2C2(2); Cor?: C,(2); Car?: 
2C2(2); Cz,?®: 2C2(2); C,,"!: 2C2(2), C2(4), 2C,(4); 
Co,!?: C,(4); Cor!®: 3C2(4); Corl: 2C2.(2), 3C,(4); 
Co": 2C2(4), C.(4); Carl’: C2(4), C.(4); Cor’: 
C2(4); Cao'®: Coo(4), C2(8), 2C.(8); Cav’: C2(8); 
Co.: 2C2n(2), 2C.(4); Cov?4: 2C2(4); Co,%: C2(4), 
C.(4); V': 8V, 12C2(2); V?: 4C2(2); V*: 2C2(2); V®: 
2C2(4); V®: 4V(2), 7C2(4); V’: 4V(4), 6C2(8); V*: 
4V(2), 6C2(4); V®: 3C2(4); Vit: 8Vi, 12C2(2), 
6C,(4); V;2: 4V(2), 2C;(4), 6C2(4); V;?: 4C>,(2), 
4V(2), 8C2(4), C.(4); Vit: 4V(2), 2C,(4), 6C2(4); 
V5: 4C>,(2), 2C2(2), 2C:2(4), 3C,(4); V2°: 2C;(4), 
2C2(4); Vi?: 4C2,(2), 3C2(4), Co(4); Va8: 2C;(4), 
3C2(4); Vi2: 4C2n(2), 2C2(4), 2C.(4); Va: 2C;(4), 
2C2(4); Vill: 2C:(4), Co(4), Ce(4); Vil?: 4C2(2), 
2C2(4), C.(4); Val: 2C2(2), 2C:(4), 2C.(4); Va": 
2C;(4), C2(4); Val®: 2C:(4); Vil®: 2C;(4), C.(4); 
Vil": 2C2n(4), Cov(4), Ci(8), C2(8), 2C.(8); Va: 
2C2n(4), Ci(8), 2C2(8), C.(8) ; Vil?: 4V2,(2), 2Con(4), 
6C2(4), C2(8), 4C,(8); Vi: 2V(4), 4C2(4), 
5C2(8), C.(8) ; Vi24: 2V(4), 4C2n(4), Cov(4), 5C2(8), 
2C,(8); Vi22: 2V(4), 2C:(8), 4C2(8); Va2*: 2V;(4), 
3Cax(8), V(8), 3C2(8), 3C2(16), 3C,(16); Vi: 
2V(8), 2C:(16), 3C2(16); V2: 4V2(2), 6C2(4), 
Ci(8), 3C.(8); Vi?®: 2V(4), 2Can(4), Ci(8), 4C2(8), 
C,.(8); V;2": 2C;(8), 3C2(8); V;,8: 4C2,(4), C2v(4), 
2C2(8), 2C.(8); Sa: 454, 3C2(2); S?: 45S,(2), 
2C2(4); Vai: 4Va, 2V(2), 2C2v(2), 5C2(4), C.(4); 
V.?: 4V(2), 2S,(2), 7C2(4); Vai: 2.S,(2), C2»(2), 
C2(4), C.(4); Vat: 2S4(2), 2C2(4); Vab: 4Va(2), 
3C2o(4), 2C2(8), 2C4(8) ; Vao: 2V(4), 2.S4(4), 5C2(8) ; 
Vai: 2S4(4), 2V(4), 4C2(8); Vai: 254(4), 2V(4), 
4C2(8); Va: 4Va(4), 2Co(8), 2C2(16), C,(16); 
Vi°: 2.54(8), 2V(8), 4C2(16); Valt: 2Va(2), V(4), 
S4(4), Cov(4), 3C2(8), C.(8); Val?: 2S4(4), 2C2(8); 
Cy: 2C4, C2(2); Ca®: 3C2(2); Ca®: Ca(2), C2(4); Ca®: 
C2(4); Cant: 4Can, 2Con(2), 2C4(2), C2(4), 2C.(4); 
Cy?: 4C>,(2), 2S,(2), 3C2(4), C,(4); Cy,*: 2S,(2), 
C,(2), 2C:(4), C2(4); Car’: 2S4(2), 2C:(4), 2C2(4); 
Cu5: 2Cu(2), Con(4), Sa(4), Ca(4), Ci(8), Cr(8), 
C.(8); CaS: 254(4), 2C:(8), C2(8); Cat: 2Cu, 
Cov(2), 3C.(4); Cav?: Ca(2), Cav(2), Co(4); Car®: 
2C2,(2), C:(4), C,(4); C,,!: C2(2), C2(4), C.(4); 
Cuy®: 2C,(2), C2(4); C..*: C,(2), C2(4); Cu’: 
3Coy(2), 2C.(4); Cas®: 2C2(4); Car®: Cav(2), Cov(4), 
2C,(8); Cal: Ce(4), Cor(4), Cr(8); Carl: Cor(4), 
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C,(8); Ca™: C2(8); D2: 4D,, 2V(2), 2C,(2), 
7C2(4); Dg: 2V(2), Cs(2), 3C2(4); De: 3C2(4); 
D4: C.(4); Da: 6V(2), 9C.(4); Df: 2V(2), 
4C2(4); D4’: 3C2(4); Dé: C.(4); Di’: 2D,(2), 
2V(4), Cs(4), 5C2(8); De®: 2V(4), 4C2(8); Da’: 
4D 4, 2V,(2), 2C(2), 7C2(4), 5C,(8); D,,”: 
2D4(2), 2Cu(2), V(4), Can(4), 2C4(4), 4C2(8), 
C,(8); DuF: 2D4(2), 2 Va(2), 2Con(4), Ca(4), Cov(4), 
4C2(8), C.(8); Dust: 2D4(2), V(4), Sa(4), Ca(4), 
C(8), 4C2(8); Da®: 2C(2), 2Vn(2), Cua(4), 
3Cov(4), 3C5(8); Dan®: 2Cu,(2), Con(4), V(4), Ca(4), 
2C2(8), C,(8) ’ D,,':2 V.(2), C4r(2), 2C2(4), C2(4), 
2C2(8), 2C.(8); D,': V(A), S4(4), C,(4), Cx(8), 
2C2(8); Da®: 4Vi(2), 2Va(2), 7C2r(4), C2(8), 
3C,(8); Da): 2Va(2), 2V2(2), V(4), Cor(4), 
4C2,(4), 3C2(8), 2Cs(8); Dax": 3V(4), Sa(4), Ci(8), 
5C2(8); Da'?: 2Va(2), 2V(4), 2Cn(4), Co(4), 
5C2(8), C.(8); Dan!®: 2Con(4), Sa(4), V(4), 3C2(8), 
C,(8); Dal: 2Vi(2), Con(4), Ss(4), 3C20(4), C2(8), 
2C,(8); Day!®: 2Va(2), 2C2.(4), C(8), Co(8), C.(8); 


Du '®: V(4), Sa(4), 2Con(4), Cov(4), 3C2(8), C.(8); 


Dy": 2D4,(2), V,(4), 74(4), C4r(4), Cx»(8), 
4C2,(8), C.(16), 3C,(16); D,,"*: D,(4), Vi(4), 
Cu(4), Vn(4), Car(8), Ca(8), 2C2(8), 2C2(16), 
2C,(16); Da!*: 2Va(4), 2Con(8), Cov(8), 2C2(16), 
C,(16); Du: V(8), S(8), Ci(16), 3C2(16); T’: 
2T, 2V(3), Cs3(4), 4C2(6); T?: 47(4), C3(16), 
2C2(24); T*: T(2), V(6), Cs(8), 2C2(12) ; T*: C3(4); 
T®: C3(8), Co(12); T,': 27%, 2Vn(3), 4C20(6), Cs(8), 
2C.(12); Ti?: T(2), 2C3:(4), V(6), Cs(8), 2C2(12); 
T)3: 27,(4), T(8), Con(24), C2.(24), Cs(32), C2(48), 
C,(48); T,!: 27(8), 2C3;:(16), C3(32), C2(48); T,°: 
T,(2), Vi(6), Csi(8), 2C2r(12), C3(16), C.(24); Ti ®: 
2C3;(4), C3(8); T;!: 2C3;:(8), C;3(16), C2(24); Te: 
2Ta, 2Va(3), Csv(4), 2C2(6), C2(12), C.(12); Ta’: 
4T4(4), Cs.(16), 2C2,(24), C,(48); Ta: Ta(2), 
Va(6), Csv(8), Sa(12), Cov(12), C2(24), Co(24); Ta’: 
T(2), V(6), 2S4(6), Cs(8), 3C2(12); Ta’: 27(8), 
25,4(24), C3(32), 2C2(48); Ta’: 2S4(12), C3(16), 
C,(24); O': 20, 2D4(3), 2C4(6), C3(8), 3C2(12); 
O2: 7(2), 2D3(4), 3V(6), C3(8), 5C2(12); O*: 
20(4), T(8), V(24), Cs(24), Cs(32), 3C2(48); O*: 
2T(8), 2D3(16), C3(32), 2C2(48); OF: O(2), D4(6), 
D;(8), V(12), Cs(12), C3(16), 3C2(24); O8: 2D;(4), 


C3(8), C2(12); O7: 2D3(4), C3(8), C2(12); O%: 
2D;3(8), 2V(12), C3(16), 3C2(24) ;O,!: 201, 2D (3), 
2Ca(6), Csv(8), 3C2(12), 3C,(24); On?: O(2), 
D4(6), Csi(8), Sa(12), Ca(12), Cs(16), 2C2(24) ; O,*: 
Ti(2), Vn(6), 2Va(6), Ds(8), 3C2.(12), C3(16), 
C2(24), C.(24); On*: Ta(2), 2Dsa(4), Va(6), Cs(8), 
V(12), Ce(12), 2C2(24), C,.(24); On®: 20,(4), 
Ta(8), Vi(24), Cav(24), Cav(32), 3C2.(48), 2C,(96) ; 
O,°: O(8), 71(8), Va(24), Can(24), C2.(48), Ca(48), 
C3(64), C2(96), C.(96); O,7: 27u(8), 2Dsa(16), 
Cs(32), Cav(48), C.(96), C2(96); O,§: T(16), 
D;(32), C3:(32), Ss(48), C3(64), 2C2(96); O,°: 
On(2), Dar(6), Dsa(8), Va(12), Cao(12), Cav(16), 
2C2»(24), C2(48), 2C,(48); O,": C3;(16), D;(16), 
V(24), S4(24), C3(32), 2C2(48); C3': 3C3; Cat: Cs; 
C3,1:2Cs;:, 2C3(2), 2C (3) ;C3:7:2Cz:, C3(2), 2C,(3); 
C;,!: 3C3v, C.(3); Ca*: Coo, C;3(2), C,(3); C;,3: 
3C3(2); Cart: 2C3(2); Cav: Cav, Co(3); Car®: Cs(2); 
D;!: 6D3, 3C3(2), 2C2(3); Ds?: 2D3, 2C3(2), 2C2(3); 
D;?: 2C2(3); Ds*: 2C2(3); Ds®: 2C2(3); Ds®: 2C2(3); 
D;’: 2D3, C3(2), 2C2(3); Daa: 2D sa, 2D3(2), Csv(2), 
2C2(3), Cs(4), 2C2(6), Cs(6); Dsa?: 3D3(2), Cs:(2), 
2C;(4), C;(6), C.(6); D3: 2D 3a, 2C3.(2), 2C2(3), 
2C2(6), C.(6); Dsa*: Da(2), Cai(2), 2Cs(4), C.(6), 
C2(6); Dsa5: 2D3a, Cav(2), 2Con(3), 2C2(6), C.(6); 
D3a°: D3(2), C3:(2), C3(4), Ci(6), C2(6); Can': 6Cs,, 
3C3(2), 2C,(3); Dat: 6Ds, 3C3(2), 2Ce20(3), 
3C,(6) ; D3)’: 3D;(2), 3C3,(2), 3C3(4), C2(6), C,(6); 
D3,*: 2D3,, 2C3.(2), Cav(2), 2C2v(3), C3(4), 3C.(6); 
D;;': D;(2), 3Cx,(2), 2C3(4), C2(6), C,{6); C,!: Cs, 
C3(2), C2(3); Cet: 2C2(3); Co®: 2C2(3) ; Co®: 2C3(2); 
Cz,!: Ces C3,(2), C2v(3), 2C,(6); Co": C;(2), C;3(4), 
C2(6); Ce: C3.(2), C3(4), C.(6); Cort: 2€3,(2), 
C.(6); Cent: 2Con, 2Csn(2), Co(2), 2C2(3), Cs(4), 
C2(6), 2C.(6); Cer?: 3Ca,(2), Csi(2), 2C3(4), Ci(6), 
C,(6); D,!: 2Ds, 2D;(2), C,(2), 2V(3), C;(4), 
5C2(6); De’: 2C.(6); De®: 2C2(6); Det: 4V(3),- 
6C2(6); D,°: 4V(3), 6C2(6); D,°: 4D;(2), 2C;3(4), 
2C2(6); De,!: 2D ex, 2D 3n(2), Cov(2), 2 Vn(3), Cav(4), 
5C2,(6),4C,(12); Don?: Do(2), Con(2), Ds(4), Can(4), 
C6(4), V(6), Con(6), Cs(8), 3C2(12), C,(12); Dg,?: 
D3,(2), Dsa(2), C(4), D;(4), C3v(4), Cx(6), 
C2(6), C3(8), C2(12), 2C,(12); Det: Dga(2), 
3D3,(2), 2C3v(4), Cx(6), C2(6), C,(12), 2C,(12). 
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The steady thermal reaction of hydrogen and oxygen at a pressure of 60 cm and a tem- 
perature in the region of 550° C was searched by absorption spectrum technique for free OH 
radicals. None were detected, although the apparatus would have been sensitive to an OH 
radical pressure of 4 10- mm. This negative result indicates that the OH radicals have a very 


short lifetime in the reacting mixture. 





I. PROBLEM 


HE existence and life history of free 

hydroxyl radicals has been investigated, as 
described in preceding papers,? by observing 
their absorption spectrum. In particular, this 
method was applied to free OH generated in 
low pressure water vapor by an electric dis- 
charge. The present paper describes an attempt 
to trace free OH in the thermal hydrogen- 
oxygen reaction. The various theories of this 
reaction agree in assuming the formation of free 
OH as an intermediary step. If we were able to 
test the concentration of free OH as a function 
of the time, interesting conclusions could be 
drawn regarding the mechanism of the reaction. 
The negative result reported below gives an 
upper limit for the free OH present under the 
most favorable conditions that can be prepared 
with a reasonable effort. 


Il. EXPERIMENT 
(a) The Thermal Reaction 


The reaction vessel, which was designed to 
have the greatest possible path length for the 
light used in the measurement of the absorption 
spectra, consisted of a Pyrex tube, of length 134 
cm, and inner diameter 4.8 cm. Heavy windows 
of fused quartz were sealed with glyptal onto 
both ends. The two gases were admitted suc- 
cessively. If the gases were admitted through an 

*At present at the Radiation Laboratory, Massa- 
chusetts Institute of Technology, Cambridge, Massa- 
chusetts. 

1 The experiments reported were performed during the 
years 1938-41. The results obtained by the first group of 
experiments were reported in Phys. Rev. 58, 1121 (1940). 


2 For references see O. Oldenberg and F. F. Rieke, J. 
Chem. Phys. 7, 485 (1939). 
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inlet attached to one end of the long cylindrical 
reaction vessel, the mixing would possibly not be 
quite thorough. In order to secure a uniform 
mixture the gases were admitted through two 
inlets properly located. 

The reaction vessel was heated by an electric 
oven 11 cm shorter than the vessel. To protect 
the seals the protruding ends of the vessel were 
cooled with water running through lead pipes. 
The temperature was between 540 and 559°C. 
Great care was taken to keep the temperature 
uniform along the vessel. This was necessary 
since the rate of the reaction strongly depends 
on the temperature,’ so that any local high tem- 
perature would cause an explosion. The temper- 
ature was kept uniform within +2°C from the 
average along a section of 80 cm length with the 
help of auxiliary heating coils wound on the same 
core with the main coil. The oven was thermo- 
stated by the arrangement described by Ban- 
croft.4 The temperature was measured with a 
platinum/platinum-rhodium thermocouple tested 
for homogeneity of the wires. To guard against 
danger from explosions the oven containing the 
reaction vessel was surrounded by copper gauze 
screens. Electrolytic hydrogen and oxygen were 
used, each being purified by passage through 
tubes containing P.O;, platinized asbestos, and 
again P.O. 

For most of the experiments, the surface of 
the reaction vessel was conditioned by successive 
reactions, separated by short periods of pumping. 
Previous work* had’ shown that this procedure 
leads to reproducible reaction rates. 

30. Oldenberg and H. S. Sommers, Jr., J. Chem. Phys. 


9, 432 (1941). 
4D. Bancroft, Rev. Sci. Inst. 13, 24, 114 (1942). 
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(b) Optical Arrangement 


The optical arrangement has been described 
in preceding papers.’ In order to make the test 
by the absorption spectrum as sensitive as 
possible, a 21-foot concave grating was used in 
the second order. The length of the spectral lines 
was reduced and so the astigmatism, which is 
caused by the concave grating, was largely re- 
moved by a cylindrical quartz lens placed in 
front of the plate. For a continuous background 
the spectrum from the positive carbon of an arc 
burning in dry nitrogen was used. This seems to 
be the most intense background available around 
3000A. For some exposures the continuous 
spectrum from a discharge through dry hydrogen 
was applied. With the apparatus available, this 
was not quite as intense as the spectrum from 
the carbon, but it was more uniform with time. 

In order to take advantage of the full resolving 
power of the grating, Eastman Process Plates or 
Spectroscopic Plates I111-O were used. The times 
of exposure were between one and one-half and 
three minutes. 

In order to increase the sensitivity of the test 
for OH and so give the negative result reported 
below greater significance, in the final group of 
experiments the light was twice reflected and so 
passed three times through the reaction vessel. 
This increased the effective length of the ab- 
sorbing gas to 240 cm. 


(c) Procedure 


The hydrogen was admitted to the reaction 


vessel first, and the oxygen second, to prevent . 


a second-limit explosion during the mixing of the 
gases. This was particularly important in the 
experiments in which the hydrogen represented 
two thirds of the mixture. Immediately after the 
admission of both gases, the reaction had its 
maximum rate; that is, the induction period was 
negligible provided that the vessel had been 
prepared by preceding runs.* The rate was 
measured by the usual method of observing the 
pressure as a function of time. In most experi- 
ments the total initial pressure of the mixture 
was 60 cm. No pressures higher than 65.4 cm 
were used, because of the danger of a third-limit 
explosion. 


5O. Oldenberg and F. F. Rieke, J. Chem. Phys. 6, 442 
(1938). 
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The absorption spectrum was photographed 
while the reaction progressed at a high rate. In 
some cases several exposures were made during 
one reaction, but in most cases only one ex- 
posure. 


Ill, RESULT 


No trace of the absorption spectrum of OH 
could be observed. The conditions were so varied 
that in every attempt a high rate was obtained 
(see preceding section). The highest rate obtain- 
able was 4.0 cm/min. In two cases it happened 
that the absorption spectra were taken at rates 
of the reaction so high that the steady reactions 
ended in violent explosions (‘‘third” explosion) 
shattering the quartz windows. Even in these 
most favorable cases the photographic plates 
showed no trace of the absorption spectrum. 
(One of these plates was underexposed, but could 
be utilized when combined with another plate 
so that the two emulsions were in contact.) 

In three experiments equal parts of hydrogen 
and oxygen were used, instead of the stoichi- 
ometric mixture, with the same negative result. 

A quantitative statement can be made on the 
basis of the recent measurements of Dwyer and 
Oldenberg,® who observed absorption spectra of 
OH with the same spectrograph and with a 
known concentration of OH. In a mixture of 
2H.O+40, at thermal equilibrium, the faintest 
observable trace of the absorption spectrum was 
found at a temperature of 768°C. From the 
equilibrium constants, the partial pressure of 
OH follows as 4.2X10-* mm. (For our estimate 
we may disregard the fact that in our case the 
conditions are a little more favorable for the 
observation of the absorption spectrum, because, 
at our lower temperature [823 instead of 1041°K ], 
the rotational energy is distributed over, not 
quite so many levels and, furthermore, the lines 
are slightly narrower.) When we take into 
account the different lengths and temperatures 
of the absorbing layers, we find that in the rapid, 
steady thermal: reaction the OH radicals, if 
present at all, are present in a mole fraction of 
less than 210-7. 

This negative result does not necessarily 
indicate that the free OH radical takes no part 


6 R. J. Dwyer and O. Oldenberg, J. Chem. Phys. 12, 351 
(1944). 
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in the thermal reaction. More plausible is the 
alternative assumption, made by von Elbe and 
Lewis,’ that, near 550°C, the rate of the reaction 
which consumes OH is so high that no appreci- 
able accumulation of OH can occur. 

At present, the theory of von Elbe and Lewis 
does not predict the concentration of OH, nor 
does any other theory, since the rates of the 
reactions producing and consuming OH are not 
sufficiently known. All schemes ever proposed 
agree in the assumption that the reaction con- 
suming OH—except adsorption at the surface— 
is OH+H:2—H.0+H. This reaction is exo- 
thermic by 15.7 kcal. Its heat of activation is not 
known. On the other hand, many different reac- 
tions producing OH have been proposed, as is 
evident in the reaction scheme of von Elbe and 
Lewis. The principal result of the present experi- 
ment is that at 550°C the reaction just men- 
tioned, consuming OH, is so rapid that no 
noticeable accumulation of OH occurs. This 
agrees with the result recently published by 
Avramenko,' who, on the basis of observations in 
dilute hydrogen flames, concluded that the reac- 
tion mentioned is very rapid. 


IV. INCIDENTAL OBSERVATIONS 
(a) Search for H,O, 


Pease® discovered that a considerable per- 
centage of H:O2 was formed at 520—550°C in 
Pyrex vessels through which a 2:1 hydrogen- 
oxygen mixture was passing at atmospheric 
pressure. The yield of HO, strongly depended 
on the treatment of the surface. Since this yield 
was as high as 1 mole H.O, to each 4 moles of 
H,0O, it seemed worth while to attempt to trace 
the concentration of H.Oz by the intensity of its 
absorption spectrum. 

It is known that the absorption spectrum of 
H,O, is continuous and has an intensity increas- 
ing toward short wave-lengths.'° Therefore it is 

02) von Elbe and B. Lewis, J. Chem. Phys. 10, 393 
(1942). 

8 L. I. Avramenko, Chem. Abs. 39, 3194 (1945); J. Phys. 
Chem. (U.S.S.R.) 18, 197-206 (1944). 

®R. N. Pease, J. Am. Chem. Soc. 52, 5106 (1930) and 
53, 3188 (1931). 

10H. C. Urey, L. H. Dawsey, and F. O. Rice, J. Am. 
Chem. Soc. 51, 1371 (1929). In this paper the absorption 
coefficients a of H2O2 are determined. However, the value 
of cd is given by a figure ten times too large. We assume 


that this is due to a misprint, so that the figures given for 
the absorption coefficients a are correct. 


convenient to use as a background the con- 
tinuous emission spectrum of the electric dis- 
charge through He. We observed the shortest 
wave-length easily photographed with a quartz 
spectrograph, that is, 2300A. The reaction 
vessel was the same as in the experiments re- 
ported above, that is, Pyrex treated by preceding 
reactions, hence presumably covered with H,O. 
The highest rate used was 2.6 cm/min. 

The result was negative. Since consecutive 
exposures were compared, the accuracy of the 
photometric method was not higher than, say, 
10 percent. This would indicate a partial pressure 
of H2O2 of less than 0.6 mm. In view of the 
results of Pease, this low limiting value should 
be attributed to the surface conditions. 


(b) Emission Spectra 


Since in the experiments described an un- 
usually deep layer of the reacting gases was 
available, it seemed worth while to investigate 
the emission spectra of the various types of the 
reaction. 

Beginning with low pressure, at the explosion 
of the second limit ten successive explosions gave 
only the faintest trace of the OH bands photo- 
graphed with a quartz spectrograph of high 
intensity. 

The steady reaction between the second and 
third limits does not emit any light of appreciable 
intensity. This agrees with the result of W. 
Jono," who needed (for photographing the OH 
bands) as many as 16,000 steady reactions at a 
rate only slightly below the first explosion limit. 

These observations are in agreement with the 
proposed reaction schemes which entail free 
radicals and atoms, but not radicals or atoms or 
molecules in excited states. 

When the first gas admitted has too low a 
pressure, inflammation (or explosion of the 
second limit) takes place at the admission of the 
other gas. When the admission is continued, 
burning takes place. The spectrum contains the 
OH bands with strongly abnormal rotation.” The 
rotational lines can be traced to about the same 
limit (quantum number 29) as in the hydrogen- 
oxygen flame. 

11 Reported by S. Horiba, Rev. Phys. Chem. Japan, 11, 


195 (1937). 
12 Q, Oldenberg, Phys. Rev. 46, 210 (1934). 
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Electrical circuits are developed to determine the normal frequencies and normal modes of 
vibration of polyatomic molecules having arbitrary configurations. The networks for each 
molecule are set up without equations from physical considerations alone and may be solved 
by an a.c. network analyzer, by numerical methods, or by digital (punch card, etc.) calculating 
machines. The networks also allow the establishment of the resultant secular equations of the 
molecule by simple inspection. The potential function considered is either a general quadratic 
form of the displacements, or it may correspond to the customary “‘valence-force’”’ system 
employing stretching, bending, and twisting of the valence bonds and valence angles. A com- 
panion paper contains network analyzer tests on the COs, OCS, and ethyl alcohol molecules 
by the “valence-force”’ model. Both the normal frequencies and normal modes of vibration 
check satisfactorily with known results within the accuracy of the instruments. 





PART I. GENERALIZATION POSTULATES 
Mathematical Models! 


O account for the fundamental frequencies 

of the infra-red and Raman spectra, the 

chemist devised several types of mathematical 
models : 

(1) The nuclei execute harmonic vibrations 
requiring a quadratic potential energy function. 

(2) The nuclei execute anharmonic vibrations 
requiring cubic and quartic, etc., terms in addi- 
tion to. the quadratic terms in the potential 
function. 

(3) The ‘‘central force’’ model assumes forces 
only along the lines connecting the nuclei. 

(4) The “valence force’’ model associates with 
the valence bonds appropriate mechanical prop- 
erties such as stretching and twisting of the 
bonds, bending of the angles between two bonds 
or between two planes containing the bonds. 

(5) The general model contains a mixture of 
the quadratic and the valence force models. 


Physical Models 


Since the analytical study of the above types 
of models becomes increasingly complex with 


* Consulting engineer. 

1 The writer wishes to acknowledge his indebtedness to 
Dr. R. M. Fuoss, Sterling Professor of Chemistry, Yale 
University for several discussions and for his encourage- 
ment to pursue the electric-circuit model study of complex 
molecules. Also he wishes to thank Drs. S. Dushman, J. P. 
Blewett, and R. P. Johnson, of the Research Laboratory 
of the General Electric Company, for their comments and 
encouragement. 


increasing numbers of atoms or radicals, attempts 
were made to introduce mechanical models.? To 
represent the valence force model, the atoms were 
replaced by weights, the bonds between them by 
springs (or beams), and the angular bending 
between two bonds was given by the bending of 
the beams themselves. Such mechanical models, 
however, do not lend themselves to satisfactory 
measurements nor do they allow the variation of 
the various weights and elastic coefficients by the 
necessary small steps. ; 

On the other hand, electric circuit models set 
up on an a.c. network analyzer can be adjusted 
to variations by rather small steps and allow 
easy measurement of both the normal fre- 
quencies and the modes of vibration not only of 
the valence force but of any type of model as 
long as it is linear. In the absence of a network 
analyzer, the electric circuits may be solved 
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Fic. 1. Masses and springs with one degree of freedom. 
(A) The OCS molecule. (B) Its electrical analog. 





2D. E. Teets and D. H. Andrews, J. Chem. Phys. 3, 
wa Murray, Deitz, and Andrews, ibid. 3, 180 
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numerically* by hand or by punch-card methods. 
The circuits also allow the setting up of the 
resultant secular equations by simple inspection. 


Systems with One Degree of: Freedom 


It has been known for over a century that a 
mass is analagous to an inductance and a spring 
to a capacitor. Let, for instance, the OCS mole- 
cule be considered, Fig. 1. The three masses 
representing the three atomic nuclei lie in a 
straight line along the x axis and are connected 
by two springs that represent the bonds between 
the nuclei. The electrical model consists of two 
capacitors connected to the ground by three 
inductors. If all the three masses vibrate along 
the x direction at one of their natural frequencies, 
in the model oscillatory electric currents of the 
same frequency flow. Then 

(1) The voltages across the inductors repre- 

sent the displacement of the corresponding 
masses in the x direction. 

(2) The currents in the inductors represent 

the forces acting on the masses. 

Two important limitations of this conventional 
electrical model should be noted. 

(1) Each mass must vibrate along one direction 
only, namely along the x axis. 

Actually in the OCS molecule the nuclei also 
vibrate at right angles to the x axis, but the elec- 
trical model is unable to represent that type of 
motion. In more general molecules each mass 
may have six degrees of freedom, namely three 
translational motions along the x, y, and z axes 
and three rotational motions around them. The 
conventional electrical model completely fails 
there. 

(2) Each valence bond also must lie along the 
x axis. 

If the angle at C is different from 180 degrees, 
the conventional electrical model again fails. Of 
course, in general molecules the bonds extend in 
every direction. 


Systems with Several Degrees of Freedom 


It is the purpose of this paper to develop elec- 
trical circuits that represent the vibration 
properties of nuclei that have six degrees of 

3 Gabriel Kron, ‘Numerical Solution of Ordinary and 


Partial Differential Equations by Means of Equivalent 
Circuits,” J. App. Phys. 16, 172-186 (March, 1945). 


freedom. They are connected by valence bonds 
that extend in space along arbitrary directions 
and are endowed with arbitrary properties. 

A mechanical system which is rather similar to 
that of vibrating molecules is a spatial frame 
structure of beams. Equivalent circuits for such 
a frame structure have already been developed 
in another publication. However, a frame struc- 
ture is not an exact analog of the chemist’s 
molecule. In some respects it is too general and 
in other respects it lacks many of the properties 
with which the molecules are endowed. For 
instance, in the molecular model there is no 
uniform bending along the length of a bond and 
any bending at the junctions takes place in a 
plane only. On the other hand, the chemist 
assumes cross-product stretching between two 
bonds, a concept that does not occur in a struc- 
tural frame. 

The equations of a simplified beam that can be 
stretched, twisted, and bent at the junctions will 
be developed here anew from fundamentals in a 
manner that will enable it to acquire properties 
that a physical beam does not possess. For that 
purpose it will be necessary to decompose each 
beam into several composite beams, each of the 
latter possessing only one particular property 
such as a stretching or an angular bending, etc. 
Each bond will then be a summation of several 
of such properties and each molecule will be a 
summation of several of such beams. 


Three Principles to Establish Models for 
Physical Systems 


The mathematical model of a physical system 
is a set of symbols arranged in the form of 
equations, while the physical model considered 
here is a set of stationary electrical coils (re- 
sistors, inductors, and capacitors) arranged in a 
configuration that spatially is analogous to the 
original system. For simple systems these models 
are also simple, but for complex configurations, 
such as a polyatomic molecule, both types of 
models soon become unmanageable and one gets 
lost in a maze of mathematical symbols or elec- 
trical coils. 

However, in going from simpler to more com- 

4Gabriel Kron, ‘“‘Tensorial Analysis and Equivalent 


Circuits of Elastic Structures,’”’ J. Frank. Inst. 238, 399- 
442 (December, 1944). 
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plex physical systems, it is possible to follow 
certain guiding principles that bound the course 
of reasoning in a narrow channel and establish 
along the flow guideposts to direct a straight-line 
progress. These principles must run a parallel 
course for either type of model and must have 
a one-to-one correspondence. They have already 
proved indispensable in establishing models for 
complex engineering systems representable either 
by sets of algebraic or by sets of partial differ- 
ential equations. 

To establish a mathematical model for a com- 
plex physical system it has been postulated that*® 


“The set of equations of a complex system composed of 
several units, each with degrees of freedom, are the 
same as the set of equations of a simple system composed 
of units, each with one (or a small number of) degrees of 
freedom except that each scalar symbol in the latter set 
is replaced by an appropriate tensor.” 


A tensor is a scalar, one-rowed matrix, rectan- 
gular or cubic matrix, in general an n-dimensional 
matrix, which obeys certain definite laws of 
transformation. A tensor is considered as the 
mathematical equivalent of a physical entity 
with nu degrees of freedom (such as a bond 
between two nuclei) and as such it is represented 
(for purposes of manipulation) by a single 
mathematical symbol. 

The parallel principle for establishing an elec- 
trical (physical) model is postulated as follows: 


“The equivalent circuit of a complex system composed 
of units with degrees of freedom is the same as the 
equivalent circuit of a simple system composed of units 
with one (or a few) degrees of freedom, except that each 
coil in the latter model is replaced by an appropriate n-wire 
transmission system.” 


The mz-wire transmission line contains in 
general as many coils as the corresponding tensor 
contains components. (See Fig. 4.) The trans- 
mission line is the electrical equivalent of a 
physical entity with » degrees of freedom and as 
such it is represented (for purposes of manipu- 
lation) by a single electrical symbol, that of a coil. 

(It is not intended to restrict the physical 
model to stationary coils and to stationary 
transmission lines, nor to restrict it to the use of 
electrical coils.) 


5 Gabriel Kron, “Equivalent Circuits of the Elastic 
Field,’’ J. App. Mech. All, 146-61 (September, 1944). 
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The connecting link between the two parallel 
postulates is a third postulate: 


“Only a tensor equation (an equation in which each 
symbol is a tensor) is representable by a physical model. 
Vice versa, all physical systems (or rather their models) can 
be mathematically represented correctly only by tensor 
equations.” 


This last postulate denies the possibility of 
constructing a physical model (an equivalent 
circuit) for any haphazard set of mathematical 
equations claiming to describe the behavior of a 
physical system. The equations must be in a 
definite tensorial form to correspond term by 
term to a model. This postulate denies the 
existence of even a mental model (the type en- 
gineers like to use, such as “‘cutting”’ of flux lines) 
dictated by a set of non-tensor equations. (The 
non-existence of a geometrical model for a non- 
tensor equation is commonly understood.) 

Since the physical system to be analyzed 
presently consists of elementary springs and 
masses whose conventional equations are already 
in a tensorial form, the third postulate will play 
only a subsidiary role in the following. (It plays 
a crucial role, however, in the model representa- 
tion of partial differential equations along curvi- 
linear coordinates or of moving and rotating 
mechanical and electromechanical systems.) The 
main emphasis will be on the first and second 
postulates to help establish comparatively simple 
equations and comparatively simple networks 
for the otherwise very complex molecular phe- 
nomena. 


The Course of Procedure 


These postulates suggest the following course 
of procedure to establish either the equations or 
the equivalent circuits of a complex physical 
system. 

(1) Divide the system into its component 
units. The collection of units will be called the 
“primitive” system. 

(2) Since the units themselves are rather 
complex, rotate and re-adjust each of the latter 
along its simplest possible reference frame (with- 
out destroying, of course, its essential unity). 

(3) Set up the equations, or the equivalent 
circuits, of each of the simplified units as if all 
the others were absent. 
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(4) Re-adjust and rotate each of the simplified 
units into the form in which it actually occurs in 
the system, either by a matrix of transformation 
(in the mathematical model) or by an actual 
rotation (in the physical model). That is, re- 
establish the primitive system of the model. 

(5) Reconnect the units (the primitive system) 
into the resultant system either by a matrix of 
transformation or by an analogous actual inter- 
connection of the equivalent circuits. 

By the use of the concept of “‘interconnection”’ 
in the form of a matrix of transformation, the 
intricate relations due to the simultaneous 
presence of several units are automatically taken 
care of without further thought. 

In really complex systems, as in organic 
molecules, the subdivisions may continue through 
several successive steps, introducing thereby 
“multiply compound”’ tensors, that is, tensors in 
which each component itself is a compound tensor 
(whose components are again tensors) and 
“‘multiply compound” networks, that is, networks 
in which each coil itself stands for a whole net- 
work each coil in which again represents, say, 
an n-wire transmission line. 


Units to be Assumed 


_ Since it is customary to assume often a mixture 
of the several types of models possible, the fol- 
lowing six types of units of a molecule will be 
studied here individually.*® 

(1) Mass of a nucleus. Occurs in all models. 
(See Fig. 2.) 

(2) Stretching of a bond between two nuclei. 
Occurs in both central and valence force models 
and it also represents the square term in a 
quadratic form (Fig. 3). 

(3) Cross-product between two bonds. Occurs 
in the quadratic potential model (Fig. 5). 

(4) Bending of the angle between two valence 
bonds in the absence of torsion. Occurs in the 
valence bond models (Fig. 7). 

(5) Torsion of a bond. 

(6) Bending of the angle between two valence 
bonds in the presence of torsion (Fig. 9). 

There are several other units possible which 
are not treated here. For instance, 

®G. Herzberg, Infrared and Raman Spectra for Poly- 


1045) Molecules (D. Van Nostrand Company, New York, 
1 , 





(7) Cross-product between two angles. 

(8) Bending of the angle between a bond and 
a plane containing several bonds. 

(9) Bending of the angle between two planes. 

When no torsion is present, each component 
of a unit has three degrees of freedom, namely, 
translations along the x, y, and z axes. In the 
presence of torsion, three additional degrees of 
freedom appear, namely rotations around the x, 
y, and z axes, denoted here by », q, and r. 


Determination of Normal Frequencies 


While the examples in the present paper 
express the equations and the electric model 
along Cartesian coordinates, the same general 
procedure applies when ,other types of coor- 
dinates are used, such as “internal” or ‘‘sym- 
metry’’ coordinates. In a more convenient refer- 
ence frame the network model may split into 
several independent networks. 

The solution of the secular equations is well 
known. The electrical network gives the desired 
information by the following steps: 

(1) Insert a constant frequency generator 
between any one of the junction points and the 
ground. 

(2) Assume an arbitrary frequency of oscil- 
lations w and adjust the value of coils (repre- 
senting the masses) accordingly. 

(3) Measure the generator current. 

As w varies, the generator current will also 
vary. Now the series of value of w at which the 
generator current becomes zero represent the normal 
frequencies of the molecule. Also, the correspond- 
ing absolute potentials at the junctions (assum- 
ing the ground at zero potential) represent the 
particular modes of vibration. 

With zero generator current the. electrical 
power stored in the network oscillates between 
the inductors and capacitors, requiring no ex- 
traneous source to maintain the oscillations. 

It may happen that for one normal frequency 
several modes are possible. In such cases placing 
the generator to a different junction produces a 
different type of mode. The measured potentials 
supply sufficient equations to determine each of 
the degenerate modes. 

It should be noted that the generator fre- 
quency is not varied as the frequency of the 
system varies. The networks are of the type 
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(A) 


Fic. 2. Equivalent circuit of a particle. (A) One-dimen- 
sional motion. (B) Three-dimensional motion. 





that require the generator frequency to remain 
constant and the inductors representing the 
masses to be adjustable. 

The measured normal frequency may be im- 
proved by Rayleigh’s principle in the following 
way: 

(1) Calculate the electrical power for the 
measured voltage distribution in all the ground 
coils, assuming w as unknown. This power is the 
kinetic energy. 

(2) Calculate the power in the rest of the 
network. This power is the elastic energy. 

The equality of the two energies determines a 
more correct w. 


PART II. THE PRIMITIVE SYSTEM 
A Particle 


The kinetic energy of a particle oscillating 
along the x axis is 7=4mx*. The equation of 
force is 

F=md*x/dt? = —w°*mx= Yx. (1) 

Expressing force as current and displacement 
as voltage, the equation of the electrical analog 
is J=YE. Its admittance is represented by a 
grounded inductor, Fig. 2(A). 

When the particle vibrates along all three 
directions, its equation is F=Y-x or I=Y-E 
where the admittance matrix is 











x y Zz 
x —w*m 
Y= -w’*m=y —wm - (2) 
Zz —w'm 

















The equivalent circuit is given in Fig. 2(B). 
In analogy with representing a tensor (or 
matrix) by a bold face letter, a compound coil 
will be represented by a coil drawn in heavy lines. 
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(A) 


Fic. 3. Equivalent circuit of a tension spring. (A) One- 
dimensional motion. (B) Three-dimensional motion (Eq. 


(7)). 


A similar matrix (along axes p, q, r) may be 
set up to represent the moment of inertia of a 
particle. It'is customary to neglect it, however, 
in molecular vibration studies. 


Stretching of a Bond 


A restoring force between two particles acting 
along the connecting line may be represented by 
a spring. Rotating the spring along the x axis 
(Fig. 3(A)), its potential energy is 


V =4k(x2—%;)*. (3) 


Differentiating with respect to x; and x2 in suc- 
cession, its equation is F=Y-x or I= Y-E where 
the admittance matrix 


Xi Xe 





Y Xi k —k (4) 
xl —k| 











Its equivalent circuit is shown in Fig. 3(A). It 
is a three-terminal (two-junction pair) network, 
or a transmission line, in which the two ground 
coils happen to be zero because of the accidental 
equality of the four k’s in Y. 

Let the spring be rotated to its original position 
whose direction cosines are /, m, and n. The 
matrix of transformation is 


1 2 








1| A, | 
Pg Bact 
2 re 





A=x l | m | n \+ (5) 
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system. (A) N-wire transmission line. (B) Its equivalent 
circuit. (C) Compound transmission line. 


The resultant equation is I’=Y’E’ where 
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1| k —k 
Y’=A,YA= |———|—, (6) 
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The resultant circuit is shown in Fig. 3(B). 
Its establishment is given in the next section. 


Equivalent Circuit of an N-Wire 
Transmission Line 


A symmetrical matrix may be looked upon, for 
the present purpose, to represent the impedance 
matrix of an n-wire transmission line, Fig. 4(B). 
(There are other possible interpretations.) The 
equivalent circuit of such a line (Fig. 4(B)) may 
be established from Y by the following steps: 

(1) The mutual admittances are represented 
by a coil between the respective terminals. A 
positive ‘‘coil” is an inductor, a negative ‘‘coil”’ 
is a capacitor. 

(2) The sum of each column, taken with op- 
posite signs, is replaced by a coil between the 
respective terminal and the ground. 

Such an n-wire transmission system is repre- 
sented symbolically by the conventional trans- 
mission line diagram drawn in heavy lines (Fig. 
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_Fic. 5. Equivalent circuit of a cross-product. (A) One- 
dimensional motion. (B) Three-dimensional motion. 
(C) Equivalent circuit of a compound coil k (Eq. (9)). 


4(C)), the “‘compound transmission line.” In 
many cases considered in this paper the ground 
coils are absent. 


The Cross-Product Terms 


A cross-product term between two bonds in 
the potential energy, such as $kAr,2Ar34, repre- 
sents a mutual spring constant between two 
springs (Fig. 5(A)) and is analogous to the mutual 
coupling between two parallel transmission lines. 

Let the two springs both be rotated along 
the x direction, still with the mutual coupling 
between them. The potential energy of the 
latter is 

V =4R(Xa2—%Xa1)(X52—X01)- (8) 


Differentiating, the admittance matrix is 


Xai Xa2 Xvi Xoo 








Xai k —k 
Xa2 —k k 
Y= ; (9) 
Xp1 k —k 
Xp2 —k k 
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The equivalent circuit of the coupling between 
two transmission lines consists of four identical 
coils k as shown in Fig. 5(A), two of them dif- 
fering only in sign. 

Let now the first spring be rotated along a line 
with direction cosines J,, mc, M2, and the other 
spring along a different line with direction 
cosines ],, ms, nm» (Fig. 5(B)). The resultant 
matrix of transformation is 











a, a: by be rv ¢ 
> a: am 
Xa1| A, | A,=X| 1, Ma Na 
| 
A Xa2 A. | (10) 
— | / / / 
Xp A, zy 
Xp2 A, | A,=x l, Mp ms 














Multiplying Eq. (9) by Eq. (10) as A:YA, the 
resultant is 


"a; as bi be 
































ai 'k —-k 
Y=" i (11) 
b,| k, —k, 
bo|—k, k, 
Xp Yo Zb 
Xa} alo lamp lane | 
k=ya| Maly Mam, Many XE 
Za| Malp NaN, NaN | 
ee ee (ee 
xX.) Yi Vo Vs 
=jJa Y, Ys Ye |- (12) 
zo Yr Ys Vo 








The transpose of k is denoted as ky. 

In the electrical network the cross-product 
term k is represented by four equal ‘‘compound” 
coils (Fig. 5(C)) just as is the single transmission 
line coupling. The compound coil k is given in 
Fig. 5(C). 


Bending of Valence Angles 


The force resisting an angular change between 
two valence bonds may be replaced, in the 
absence of torsion in the molecule, by an angular 
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Fic. 6. Equivalent circuit of an angular spring between 
two distant bars. (A) One-dimensional mofion. (B) Three- 
dimensional motion. 


spring connecting two bars having no rigidity. 
Let the two bars, @ and 6, with the angular 
spring, be rotated along the x axis (Fig. 6(A)). 
The angle by which one of the bars is displaced 
is (Yao — a1) /La; hence the potential energy stored 





(C) 


Fic. 7. Equivalent circuit of an angular spring between 
two neighboring bars. (A) Angular spring. (B) Compound 
network. (C) Actual network (Table I). 
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in the spring is 
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The equivalent circuit of the angular spring 
is two-coupled transmission lines, Fig. 6(A). The 











































































































1 [¥a2—Yar Yoo—Yu f : 
=> ieee Pape (13) coupling has a sign opposite to that in Fig. 
Fi 5(A). 
Differentiating, the admittance matrix is In place of the directions of the bars themselves, 
Yai Yar Yo Yor let the directions of the lines perpendicular to the 
~ a ee Pe ae bars and lying in their plane be considered and 
let their direction cosines be denoted by J,, ma, 
ware —he Ra Rab | —Rar (14) "a and J, ms, m», respectively. The resultant 
Yui| — Ras Rab ky | —ke transformation matrix A is the same as in Eq. 
Fal ha| —ka| —hs ks (10) except on the left-hand side the unit vector 
x is replaced by y. The resultant Y’ is by 
where kg=k/L.2; kp=k/Le?; and kas=k/LaLls. AYA. 
ai a2 b; be 
a) ok, | —ke | —kov | kw | 9% P 
a| —k&, | kK | Ke | —Kes | alY. | Yas! 
Y= = Saat (15) 
b; — Kav: Kav: k, —k, b » = Y, | 
be Kav: — Kas: —k, k, | 
x y 2 junction is 
x lale L.my Ln» ay’ ay’ b.’ 
Kin=Y| Maly | MaMy | Mano LL, (16) ai 1 
Z Nal NaMps NaNdb ao 1 ~ 
a" ' , (17) 
k, (and k,) are the same as kj» with all “db” sub- . , 
scripts replaced by ‘‘a’’ (or ‘“‘a” replaced by “‘b’’). be 1 
The compound electrical network is given in 
Fig. 6(B). It is the same as Fig. 4(B) except that > 2: 
the mutual coils k,, have opposite signs. xi 1 
Angle between Two Neighboring Bonds l=y 1 (18) 
When the two bonds forming a valence angle J 
are connected together (Fig. 7(A)), then one of 
the junctions disappears. The connection matrix By the law of transformation Y’=A,YA, the 
connecting two junctions a; and }; into one resultant admittance matrix is 
ay’ ay’ b,’ 
ay| Ket+k,—Kas—Kase  —KatKuas: —k, +k, 
Y” =a,’ —k,+K., k, — ka», (19) 
b,’ —k,+Kas: — Kase k, 











The equivalent circuit is given in Fig. 7 and in Table I. In two dimensions, ng=n,=0. 
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Torsion of a Bond 


The treatment of the torsion of a bond is 
analogous to that of stretching except that the 
variables and unit vectors x, y, z are replaced 
by p, q, r. The potential energy is V= $k(p2—p1)’. 
When the bond lies along x, its Y is given in Eq. 
(4) and when lying in a general direction, its Y 
is given in Eq. (6). The equivalent circuits are 
similar to Fig. 3. 

It should be noted that the previously de- 
veloped circuits cannot be interconnected: with 
the circuit representing torsion, or rather there 
would be no interconnection between them. The 
bending of valence angles in the presence of 
torsion, newly treated in the next section, allows 
the interconnection of the x, y, z and p, q, r axes. 


Bending and Torsion of Valence Angles 


In the presence of torsion the single angular 
spring with a force constant k connecting two 


“TABLE I. Admittances of coils in Fig. 7(C). 
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Y,=k Y;=k——_ Y33= Yi9= 
L.* Ly? LeaLy Lal 
Ma? my? Moms Mand 
Y2,= Ys=k Yyu= Yoo=k 
LS L;? Lela LoL 
Ne? ny? NaN MbNa 
Y;=k Yo=k Yi;=k Yu =k 
e L;* Lela LaLy 
laMa lymy lamp 
Y,=k Yio= Yis=k _ 
a? L;? LoL 
latte lyns lyma 
Y;= Yu=k Y= 
LL." L;? Lela 
MaNa MyNns lan 
Y,=k- Y= Yis=k A 
a? L;? Leks 





Y= Y4+ Yio— Yis— Yur Ya=— Yet+ Yu 
Yo=V¥stYu—Yis—VYio9 Vi=— Yut Vio Y,= — Y3+ Vis 


Yo=—Yst+ Vis 


Ye= Vet Via—VYis— Yio YVi=—Vst+ Yuu Ye=— Yiot Vir 
Ye=—Vit Nis Ye=—YVutVa Yr=— Yiot Vis 
Ye=— Vt Yur Yi=—¥2t+Vis Ys=— Viet Yoo 
Y;=— Vit Vic Ym=—YstVis Yi=—Yut Vis 
Yo=— ¥ot+ Yuu Yn=—YVetYoo Yu=—Vot Vis 








(A) (B) 


Fic. 8. Replacement of angular spring by two beams with 
torsion. (A) Angular spring. (B) Two beams. 


bars has to be replaced by two springs in series, 
each with a force constant 2k, Fig. 8. Each spring 
now resists the angular change of a single bar so 
that each bar can be treated independently of 
the other. , 

It should be noted that when several angles 
originate on a bond, the bond is considered as 
composed of several bars in parallel, each bar 
being associated with a different angle and a 
different angular force constant k. (There may 
also be two more bars belonging to the same 
bond, relating to a stretching and a twisting 
force constant.) 

Let the bond be rotated along the x axis (Fig. 
9(A)). The angle by which the bar is displaced 
is (ye—1)/L—n, where 1; is the angular rotation 
of end 1 around the z axis. The potential energy 
in the spring is 











V=4(2k) (“—*-1). (20) 
iL 
Differentiating, the admittance matrix is 
yi y2 YT 
2k 2k 2k 
"> pl T 
y 2k 2k 2k (21) 
"| Bo L ' 
2k 2k 
Ti — an 2k 
} a L 














Its equivalent circuit, Fig. 9(A), now contains 
ground coils. 

Let the bond be rotated into its actual position 
and let the following directions be defined : 

(1) The direction cosine of the line per- 
pendicular to the bond and lying in the plane of 
the two bonds (the plane of the bending) be 
lh, mM, NN}. 











28 GABRIEL KRON 














1W—— 2 
-2k/L? 
”; it 
2k/L tS 
1 -2k/L 
-2k/L = = -2k -2k/L 
je 
(A) 


























(B) 


Fic. 9. Beam equivalent of an angular spring in the presence of torsion. 
(A) One-dimensional motion. (B) Three-dimensional motion. 
















































































(2) The direction cosine of the line perpen- P q J 
dicular to the plane of the two bonds be Lo, ma, N2. . hls heats lem 
The transformation matrix is 
ky=y) mle mM m\Nz I" (24) 
zs ¥ ss 
Z s &, poe Z| Nylo N\Moe nn. | 
yi| Ai | wees  ™ ; | ; — , 
a A | k, (and kz) are the same as ky, with all “2” sub- 
pe Rist = Ba & ft (22) scripts replaced by “1” (or “1” replaced by ‘‘2”). 
| A» ae | Its equivalent circuit is shown in Figs. 9 (B) and 
| A> =T1T)| ls M2 Ne | P 
—, = | 10 and in Tables II and III. 
The resultant ¥ is "PART III. THE RESULTANT SYSTEM 
1 2 1, The Resultant Physical Model 
2k, 2k, 2ki> Let the four-atom molecule of Fig. 11(A) be 
os ar L considered. Its primitive system is shown in Fig. 
i z ‘ 
| | 11(B). There exist a cross-product term between 
- 2K 2k 2K» the two distant bonds f and 7, and a valence 
Y=2} — - , (23) : : 
Lc? 2 i angle a between two neighboring bonds. The 
2k: 20, | potential energy of the system is 
1, L my L 2K2» | V=3[Rav(Aras)?+Roe( Aree)? +Raa(Araa)* 
— | +k(Araa) (Arse) +Ra(Aa)?]. (25) 
2 i * Fh if 22 
' SS ¢ My 
aaa es = k | 
: Zz OP )/ | 
x) 4} Xo 
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Fic. 10. Equivalent circuit of a beam. (Tables III and IV.) End 1 rigid, end 2 
free to rotate. 
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TABLE II. Admittances of coils on Fig. 11 lying between various junctions. 




















wy = Vs | Y2>22 = Yo Xe—>y2= Vo 1X2 = — Vox 
a=Y; yz = Ye | a= VY, yo= — Vu 
pi= Vio pi= V2 21> Pi = Vee | zo= — Vo 
gi = Yx n= Yeo g= Yo Pir = Vie | 
r= Vos ri= Vu r= Vu n=Vi7 nan = Vie 
X2=— VY; X2=— Vy x2=— Ys x2= — Vig x2= — Yo 
y= — Ve n=— Vs yw=— Ve y= — Vos ¥2= — Yoo 
z2=—YV; z2=— Ve z=— V3 22= — Vo zo — Vor 





ground—>x; = — (Yio-+ Yoo+ Yes) 
y= —(Yos+ Yoot You) 
a= — (Yeo+ Yer + Yuu) 
p= —(¥r+ Viet Yu) 
qi = —(Yiot ¥st+ V2) 
r= —(VYut Vit Yo) 


ground—>x3= Vio+ Yo2+ Yes 
2 = Yos+ Yoot+ Yo 
y= Yost Yor+ Yo 











Each expression in V represents a unit in the 
primitive system. 

The resultant electrical model, Fig. 12(B), is 
built up from the compound coils, Fig. 12(A), in 
a manner analogous to the interconnection of 
the units of the original system. The step from 
the compound coils to the actual coils consists of 





interconnecting all analogous x, y, z or p, q, r 
terminals separately as shown in Fig. 13. 
The physical interpretation of the network is 
as follows: 
(1) Relating to each particle 
(a) The voltage from a terminal to 
ground gives the total displacement 


TABLE III. Definition of admittances in Table II. 
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(A) (B) 


Fic. 11. Assumed molecular system X YZW. (A) Resultant 
system. (B) Primitive system. 


(linear or angular) of the particle in 
the corresponding direction. 

(6) The current between a terminal and 
a ground gives the resultant force or 
torque acting on the particle in the 
corresponding direction 

(2) Relating to each valence bond 

(a) The difference of potential between 
two terminals gives the change in the 
length or angular position of the bond 
in the corresponding direction. 

(b) The current between two terminals 
gives the corresponding force in the 
bond. 

(3) Relating to the total system 

(a) The total power in the ground coils 
representing the masses gives the total 
kinetic energy of the vibrating system. 





KRON 














(B) 


Fic. 12. Electric circuit model. (A) Primitive system. 
(B) Resultant system. 


(6) The total power in the rest of the 
network gives the total potential 
energy stored in the vibrating system. 


The Resultant Mathematical Model 


The units of the mathematical model consist 
of the Y’ matrices of the four masses (Eq. (2)), 
three stretching springs (Eq. (6)), one cross- 
product term (Eq. (11)), and one angular spring 
(Eq. (15)), all rotated in their proper spatial 
positions. These matrices are arranged along the 
diagonal line of a compound Y’ (Eq. (26)) with 
thirteen rows and columns representing the 
equations of the primitive system I= YE, namely 
thirty-nine scalar equations. ° 

The interconnected system of Fig. 12(B) has 
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a| —w’m, | | 
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only five junctions, hence four-junction pairs. If 
each of the néw junction pairs is selected so as to 
include the ground, the relation between the 
voltages (displacements) of the primitive and 
the resultant system is the matrix of trans- 
formation A, Eq. (27). 























— a 


(A) (B) 


























/ / , / 
a’ bc d Fic. 13 Interconnection of coils. (A) Interconnection of 
ly l ] compounds coils. (B) Interconnection of actual coils. 
as | 

| 
b| I where I is the unit matrix of Eq. (18). The ad- 
P 2 -_ mittance matrix of the resultant system is by 

Ce ee ae A/Y’A’. 

d) I The equations I’’= Y’’E”’. represent the twelve 
e:| I | | secular equations of the molecule. 

= 2-5 ae oe 
a, Ty | Relations between the Two Models 

A’=f | I | |, (27) . 

y ee Fe Te Bee The secular equations I= YE may be estab- 
f | |i _ lished by a simple inspection of the electric 
gil; | | | circuit model of Fig. 12(B). In particular, Y”’ is 

| Tr = established as follows. 

2 | | . . 

Lieder (1) The mutual admittance terms are the coils 
hy} | _ {ft connecting corresponding terminals. 
ij I | | (2) The main diagonal terms are the sum of 
i |/— = TF the admittances entering each junction, taken 

Ff fj. jy with opposite sign. 

a c d 
,' an — 
| l 
al YetYatYi—w'm, | Ysie-Yet+¥ou =| —Ysie | = Yos—Yi 
Y”’ b| Yei— Yee t+ Yooae | Y.+Y¥;+Y,2—w*m, —Yy, | —Ysit Yoo (28 
= 25) 
c —Y;; —Yy Y;—w’m, Yy; 
d Yost — Y; —YVyittYoo: Yyit seadicllihdiosens 











Vice versa, the electric circuit model may be 
established from the secular equation by the 
reverse steps of those given previously. The two 





methods of {setting up the electrical model may 
serve as a check on the correctness of the net- 
work. 
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Network Analyzer Tests of Equivalent Circuits of Vibrating 
Polyatomic Molecules 


G. K. CARTER* AND G. Kron** 
General Electric Company, Schenectady, New York 
(Received August 20, 1945) 





To illustrate the use of the equivalent circuits derived in the companion paper for vibrating 
polyatomic molecules, this paper describes network-analyzer tests of circuits for three suc- 
cessively more complicated molecules. The first is the symmetrical linear CO, molecule; the 
second is the unsymmetrical linear OCS molecule; and the third is the unsymmetrical non- 
linear CH;CH2OH molecule. In each case, starting from given constants of the molecular 
structure, the normal frequencies and modes are measured by means of the equivalent circuit. 
It is shown that for these cases the analyzer measurements agree satisfactorily with the known 


normal-frequency characteristics. 





NETWORK ANALYZER 


HE alternating-current network analyzer! 

consists of a set of adjustable resistance, 
inductance, and capacitance units, each con- 
nected to a pair of flexible cords and plugs. Con- 
nections between units to form any desired 
network are made by inserting the plugs in 
adjacent jacks in a jack panel. As many units as 
desired can be connected to a common point. 
Alternating-current electric power is supplied by 
a motor-generator set to individual generator 
units so that several different voltages, inde- 
pendently adjustable, can be inserted into dif- 
ferent parts of the network. A centrally located 
set of measuring instruments can be connected 
to any unit or circuit by means of a set of key 
switches. 
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Fic. 1. Equivalent circuit of the CO» molecule. 











* Analytical Division, Central Station Engineering 
— General Electric Company, Schenectady, New 

ork. 

** Consulting Engineer, General Electric Company, 
Schenectady, New York. 

1H. P. Kuehni and R. G. Lorraine, ‘‘A New A-C Net- 
work Analyzer,” Trans. A.I.E.E. 57, 67 (1938). 
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Procedure 


The equivalent circuits for the CO, OCS, and 
CH3;CH:OH molecules are shown in Figs. 1-3, 
respectively. These circuits were obtained by 
methods described in the companion paper,’ using 
the molecular constants listed in Table I; the 
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Fic. 2. Equivalent circuit of the OCS molecule. 








Fig. 3. Equivalent circuit of the ehtyl-alcohol molecule. 


2G. Kron, “Electric Circuit Models for the Vibration 
Spectrum of Polyatomic Molecules,” J. Chem. Phys. 14, 
19 (1945). 
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Fic. 4. Variation of generator current with frequency, 


constants are based on previous data,*4 with the 
angular spring constants of the former reference 
corrected by a factor of one-fourth. 

The process of locating the normal frequencies 
from these circuits is similar to that already 
described in connection with the study of the 
Schrédinger equation.® The circuits specify both 


TABLE I. Molecular constants. 








CO, L=1.16 k=15.3 m-=12 
ka= 0.772 Mo= 16 

OCS Lee 1.16 Roe= 14.2 m, = 12 
Le 1.54 ka= 8 Mo = 16 

kaz= 0.655 m,= 32 

CH;CH2OH L,=1.37 k= 4.1 Moy; = 15 
L2=1.54 ko= 3.6 Moy, = 14 

a= 109.5° ka = 1.102 Moy = 17 








The stretching force-constants k are expressed in 105 dyne/cm, the 
bending force-constants kg in 10-" dyne-cm/radian; m is the atomic 
weight (O =16); and the wave-length »v =1303w. 


3 W. G. Penney and G. B. B. M. Sutherland, ‘On the 
Relation Between the Form, Force Constants, and 
Vibration Frequencies of Triatomic Systems,” Proc. Roy. 
Soc. 156, 654-678 (1936). 

* Paul C. Cross and J. H. Van Vleck, ‘“‘Molecular Vibra- 
tions of Three Particle Systems with Special Applications 
to the Ethyl Halides and Ethyl Alcohol,” J. Chem. Phys. 
1, 350 (1933). 

5G. Kron and G. K. Carter, ‘‘A-C Network Analyzer 


positive and negative admittances; since nega- 
tive-resistance units are difficult to work with, 
positive and negative reactive admittances are 
used, with a.c. excitation. The circuit is excited 
by applying a generator at any junction point, 
and the forced-vibration frequency is varied by 
varying the admittances of the units expressed in 
terms of w. For any given w, the generator current 
represents the applied force required to force 
vibration at that frequency; so that when this 
current is plotted as a function of w, as in Fig. 4, 
the points at which the resulting curve passes 
through zero should correspond to frequencies at 
which no external force is required to sustain the 


TABLE II. Measured and calculated normal 








frequencies. 
Vi v2 Vs 
CO» Test 1273 670 2447 
theory 1274 668 2440 
OCS Test 860 526 2072 
theory 859 527 2079 
CH;CH,OH Test 897 1009 421 
theory 910 1020 430 








Study of the Schrédinger Equation,’’ Phys. Rev. 67, 44- 
49 (1945). 
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67.5 





100 


























-84.5 


Fic. 5. Measured normal modes of vibration of 
the ethyl-alcohol. 


vibration, that is, to the natural frequencies of 
the circuit and of the molecule which it repre- 
sents. Ideally, the curve shown should be a con- 


tinuously increasing function with increasing w, 
passing through a series of infinite discontinuities ; 
these discontinuities represent frequencies for 
which the point of excitation is a nodal point; 
practically, of course, the current supplied is 
incapable of increasing without limit, and instead 
it changes very rapidly from positive through 
zero to negative before once more taking the 
upward course. These extra zero points are 
generally obvious; and in any case the true 
natural frequencies can and should be checked 
by using the exciting generator at more than one 
junction of the network. 

For the three molecules being considered, 
Table II shows how the natural frequencies 
measured on the network analyzer compare with 
those calculated from the constants in Table | 
(omitting the obvious case of zero frequency). It 
is seen that the check is very good; in fact, the 
discrepancies generally are less than might have 
been expected from uncorrected readings with 
the equipment used. 

Not only the natural frequencies, but also the 
natural modes, can be read from the circuits by 
means of the analyzer. The latter are given 
simply by the voltages of the several junction 
points at the natural frequencies. Figure 5 shows, 
for example, the natural modes of the CH; CH2OH 
molecule as measured in this way; these are 
similar to the corresponding modes shown in 
reference 4. 





peeeey Yn ee ne 


————— 


wi ae aes 


i 
le 


errr per: 
aa SS 9 SRR Vet SR PSEA 


nee 


te Naam 


THE JOURNAL OF CHEMICAL PHYSICS 


VOLUME 14, NUMBER 1 


JANUARY, 1946 


The Asymmetric Rotor 


III. Punched-Card Methods of Constructing Band Spectra 


GILBERT W. KING, Arthur D. Little, Inc., Cambridge, Massachusetts, 
Paut C. Cross, Brown University, Providence, Rhode Island, 


AND 


GrorGeE B. Tuomas, Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received October 10, 1945) 


The stochastic method of analysis of band spectra of asymmetric rotors involves the calcula- 
tion of line strengths and positions for successive estimates of six molecular constants, until the 
calculated intensities at all wave-lengths agree satisfactorily with the observed. This paper 
shows that these extensive calculations, as based on tables previously published, are well 
adapted to standard punched-card equipment. In‘addition to the above calculations, it is 
possible to construct the appearance of the rotational structure for various values of slit 
width. The procedure is carried through for HS with some approximate constants. Very satis- 
factory representations of observed spectra are obtained for completely resolved spectra, as in 
the photographic region; for slit width 1 and 2 cm™, as in the near infra-red; for slit width 


10 cm, as in poorly resolved infra-red spectra. 





I, INTRODUCTION 


AND spectra still remain an important 

source of information about molecular 
structure and of precise values of the constants 
for calculating the thermodynamic properties of 
simple substances. The analysis of bands of 
polyatomic molecules is very difficult when the 
three moments of inertia are different, that is 
when the molecule is an asymmetric rotor. 
Figure 1 shows typical infra-red bands, in this 
case of H.S.' Although these bands are not com- 
pletely resolved into lines, one feels that there is 
enough information given by the position and 
intensities of the peaks to find out the dimensions 
of the molecule. However, only a very few of 
such bands have been analyzed. The difficulty 
arises partly because the positions of the peaks 
in the case of the asymmetric rotor are very 
irregular. It seems that the only approach to 
interpretation of these bands, whether more or 
less resolved than the ones shown, is a stochastic 
one. One assumes the dimensions of the molecule 
in the upper and lower states, calculates the 
position and intensity of the lines, and compares 
with the observed. This has been done success- 


1H. H. Nielsen and E. F. Barker, Phys. Rev. 37, 727 
(1931). A reproduction of the completely resolved band of 
H2S at 10,100A was given by P. C. Cross, J. Chem. Phys. 
5, 370 (1937). 


fully for only a few completely resolved bands 
of asymmetric rotors, namely, a number of bands 
of H2O and one band of H2S, and even in these 
cases the comparison of intensities has been 
based on symmetric-rotor approximations. 

The position of a line in a band is given by 


Vor att; I’'.e= yt+E'(a’, b’, Chee 
—E"(a", b’’, cl") gu, ry (1) 


where vo is the band center, E(a, b, c)s,, being 
the rotational energy of the state defined by 
J,7; double primes indicate the lower state, 
single primes the upper state. In the case of the 
asymmetric rotor the rotational energies can be 
determined only by solving polynomials of 
various high degrees. In a previous paper’? we 
believe we have simplified this part of the 
problem a great deal by tabulating a quantity 
E(x) in terms of which the rotational energy, a 
function of the three reciprocal moments a2 b2 c, 
can be expressed as 


E(a, b, c)s,-=3(a—c)E(k)ys,- 
+3(a+c)J(J+1), (2) 


where «=(2b—a—c)/(a—c) is a parameter of 


asymmetry varying from —1 to +1, and E(x),,,, 
called reduced energy levels, have been com- 


2G. W. King, R. M. Hainer, and P. C. Cross, J. Chem. 
Phys. 11, 27-42 (1943), hereafter referred to as I. 
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Fic. 1. Absorption bands in hydrogen sulfide. A, Prism curve with low dis- 
persion. B, Fine structure in the band at 2.6u. C, Fine structure in the band at 


3.7. (Copied from reference 1.) 


puted* up to J=12 by solving continued frac- 
tions for all possible values of « in intervals of 0.1. 
Thus the energy levels for any molecule can be 
computed by (2) from the table of E(x) by 
simple multiplication and addition. 

The intensity of a line in an absorption band 
as obtained‘ from I]—(1) is 


Tgerget; J’, =const v(1 —EmVEP) 8 ary ar 
Xexp—E"(a"’, b,c") y,¢/RT 
ne ze Do Brg | 2a, 2, M’'": J’,7',M’, (3) 


F=X,Y,Z M’’ M’ 
where » is the line position (1) in wave numbers, 
k being the Boltzmann factor in cm~ deg, and 
gu,- is the nuclear-spin weight factor. The 
constant factor is given explicitly by II—(1) 
and I]—(2). It includes the partition function 
and the absolute magnitude of the permanent 
and induced electric moment. The last term is 
the square of a matrix element of the direction 
cosine between the radiation field and molecular 
axes. The subscript g refers to the direction of 
the permanent and induced components of the 
electric moment in fixed molecular coordinates. 
The triple sum is the strength of a line whose 
components are unresolved in an isotropic radi- 
ation field with no external field. In II we have 
computed the values of the line strengths for 


3 Values of E(x), x=—1.0 (0.1) 0.0, for J=11 and 12 
will be published shortly. 

4p, é& Cross, R. M. Hainer, and G. W. King, J. Chem. 
Phys. 12, 210 (1944), hereafter referred to as II. 


five values of x(+1, +0.5, 0), which adequately 
cover the whole range of asymmetry possible. 

Equations (1) and (3) are to be evaluated for 
all important lines of the spectrum. Even with 
the tables of energy levels and line strengths, the 
construction of a spectrum, which is composed 
even in the simplest cases of several hundred 
lines, is a considerable undertaking. The com- 
putational labor can be reduced by -making use 
of certain calculating machines. 


II. ACCOUNTING-MACHINE METHODS 


It is the purpose of this paper to illustrate the 
use of business accounting machines in carrying 
out certain types of quantum-mechanical calcu- 
lations, in this case the construction of spectra by 
the evaluation of the above equations for line 
position and line intensity. Punched-card ma- 
chines are probably the most suitable for this 
kind of work, for although they are only capable 
of adding, subtracting, and multiplying, they are 
well suited to the quantum-mechanical opera- 
tions such as symmetry classifications, sorting 
according to parity or changes in quantum 
numbers, etc. 

In the International Business Machines all the 
quantities used in the calculations are supplied 
on cards, the numbers or identification letters 
appearing as holes punched in certain locations. 
The machines read the information by brushes 
which make electrical contacts through the 
holes as the card runs through the machine. 
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Several types of machines are available whose 
principal operations are briefly described below 
with special reference to the construction of 
band spectra.® 

The ultimate object of this work is to represent 
a spectrum by a set of cards, each card corre- 
sponding to a line and bearing the line position, 
line intensity, the quantum numbers of the two 
states, and the symmetry classification. These 
cards can be prepared and handled by the fol- 
lowing machines. 

Key Punch—The original data can be punched 
onto cards by hand from printed tables. Since 
this step is subject to errors of the operator, the 
amount of material punched by hand should be 
kept to a minimum. It is necessary that the 
information on the cards be thoroughly checked. 
Various machine methods of testing the accuracy 
are described below. 

Sorter—Fundamentally this machine sepa- 
rates a set of cards into ten groups depending on 
the value of a digit in any one column. It can be 
used to arrange the cards representing transi- 
tions in order of increasing wave number or to 
pick out the lines of intensity above any chosen 
value. It is also possible to classify the lines, for 
instance, into sub-branches, by sorting on the 
changes in quantum numbers. 

In addition to the field in which digits are 
punched there are X and Y zones which can be 
used in combination with the digits to indicate 
letters,-or used alone to indicate sign, parity, 
etc., or for control purposes. Thus lines can be 
sorted into symmetry classes determined by the 
parity of certain quantum numbers. 

Collator—This machine is used to combine 
two sets of cards so that a card of the “primary”’ 
set defined say by a certain quantum number 
precedes all the cards in the ‘‘secondary”’ set with 
the same quantum number. 

Reproducer—When two sets have been col- 
lated, the combined set can be put through this 
machine so that information is reproduced from 
the leading card (from the primary set) into all 
the cards with the samé quantum numbers of the 


5A description of punched-card machines and of their 
application to astronomial calculations can be found in 
W. J. Eckert, Punched-Card Methods in Scientific Computa- 
tion (Thomas J. Watson Astronomical Computing Bureau, 
Columbia University, New York, 1940). 


secondary set. This operation is called “gang 
punching.”’ 

Multiplie-—This is the machine on which 
addition, subtraction, and multiplication can be 
done. The two numbers which are to be combined 
must appear on the same card, the answer being 
punched elsewhere on the card. Combined 
operations, such as AXB+C, can be done 
simultaneously in this machine, the addition 
being called ‘‘cross-footing.”’ 

Group multiplying is the operation of multi- 
plying all cards put through the machine by the 
same quantity which is inserted in the machine 
by a leading specially punched card. It is then 
not necessary to use up any columns of the 
detail card for the multiplier. 

It is desirable in scientific work to have the 
Multiplier modified to take into account different 
signs of the terms in an equation such as (1); 
otherwise, the cards have to be sorted according 
to the signs and handled separately: 

Summary Tabulator—At any stage in the 
operations the cards can be put through a 
Tabulator wired so that any or all of the in- 
formation can be printed. Thus, the final results 
can be listed in any convenient form. 

In addition, the Tabulator can sum the quan- 
tities on any part of the card and print either the 
grand total or subsidiary totals, the printing of 
which can be controlled by changes in identi- 
fication columns, e.g., when a quantum number 
changes or when the tens digits change. Sum- 
mation can also be controlled by X zone punches. 


Operation 


In conclusion, it should be pointed out that 
all the steps in the calculations can be set up in 
a standardized form, and the actual preparation 
and handling of the cards is done by specially 
trained operators at some computing center. All 
the spectroscopist has to do is to supply the 
necessary eight constants, namely, three recip- 
rocal moments of inertia for both the upper and 
lower states, the absolute temperature, and the 
band center. 


Punched-Card Table of Reduced 
Energy Levels 


It was not convenient to obtain the reduced 
energy levels E(x) by punched-card methods 
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since these were calculated by continued frac- 
tions, and punched-card machines are not 
adapted to frequent division. The table of 
reduced energy levels (I—Appendix) was set up 
on punched cards on the Key Punch. One card 
was made for each level for each value of x 
(running from —1 to 0 in intervals of 0.1), bear- 
ing the values of J, K_1, Ki, and E(x). 

The 1-values of the levels were then obtained 
by finding the difference s=K_,—K;, on the 
Multiplier. The cards were sorted according to 
the parity of the K’s and the resulting four 
symmetry species identified by punching the 
characters ee, eo, oe, 00 (I—Sec. IV). Similarly 
by sorting on the last digit of tr and X punching 
the odd numbers, the parity of + was indicated 
on the cards for use in applying weight factors. 
The parity of y was obtained from that of J+7 
for use in grouping into sub-branches. 

The major source of error in punched-card 
operations is the copying of tabulated material 
onto punched cards, which is done by an oper- 
ator. The punched cards representing the table 
had to be carefully checked. This was done very 
easily as follows: the cards were sorted on J and 
the four symmetry classifications and then put 
through a tabulator which printed the sums of 
each set with the same J and symmetry. These 
numbers are the traces of the sub-matrices whose 
values are known.® 


Calculation of Term Values 


To obtain the actual energy levels of any 
molecule, the required information is copied from 
the master set described above. In general it is 
necessary to interpolate to find E(x) for inter- 
mediate values of x. This can be done by machine 
methods since interpolation involves only multi- 
plication and addition. It is sometimes more 
convenient, however, not to interpolate, but to 
carry through the calculation for line position 
v for several values of «x, and interpolate the line 
position. In many cases it will be found that the 
variation of » over the interval of 0.1 in «x is 
either negligible or of the order of the centrifugal- 
distortion correction. This method emphasizes 
the sensitivity of certain lines to the asymmetry 
parameter (which in XY». type triatomic mole- 


6 [—Table X. 
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cules is a function of the valence angle only, not 
of the interatomic distances). The positions of 
these lines can be used to determine «x accurately. 
The original table of reduced energy levels 
(I—Appendix) covered only the range of nega- 
tive x. Values for positive x are obtained by 
means of the symmetry relation (I—12): 


E(«)s,x_,.K,ir= —E(—«)s,K,K_,:-2- (4) 


The inversion of K_; and K, and change of signs 
of 7 can be done automatically on reproducing 
from the master set by interchanging appropriate 
wires. 

The value of $(a+c)/J(J+1) is punched onto 
the appropriate cards by collating and gang 
punching, a process to be described in more 
detail in another connection. The pack of cards 
is put through the Multiplier, in which the 
factor }(a—c) is introduced as a group multiplier. 
The cards run through the machine so that the 
product $(a—c)E(x) is obtained and simul- 
taneously added to $(a+c)/J(J+1), so that (2) 
is evaluated at once and the answer E(a, b, c) 
punched in other columns of the card. 

This whole operation is repeated for the chosen 
constants of the upper state. 


Punched-Card Table of Line Strengths 


Another master set of cards has been made, 
one card for each possible transition between 
rotational levels, J <13, for all three components 
of the electric moment, bearing the quantum 
numbers, J’’, K’’_; and K”’; of the lower state, 
and J’, K’_; and K’; of the upper state, together 
with the line strengths. 

Since the line strengths were obtained by 
ordinary methods of calculation, they had to be 
key punched onto cards by a hand operation. 
This introduced numerous errors of identification 
which were fairly easily corrected by the appli- 
cation of sum rules (II—20 to 24). The cards 
were sorted by direction of moment, level iden- 
tification’ J’, K’’_,, K”,; (or by J’, K’_4, K's), 
and the sums of the line strengths obtained on 
the Tabulator. These sums. being small whole 





7 Since the double suffix K K-notation involves 4 integers 
(J <13) and the r-notation involves 2 integers plus a sign, 
one sorting operation can be eliminated by the use of J+r 
(running from 0 to 2/+1) as the index of the levels with 
the same J. 
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numbers, errors could be detected by visual 
inspection. The use of three different sum rules 
finally eliminated all errors except those due to 
rounding off. 


Punched-Card Table Representing a Band 


A copy of this master set of line strengths was 
made for the chosen value of « and direction of 
the electric moment. Interpolation may be 
necessary, although, as in the case of E(x), it is 
sometimes of interest to carry through the cal- 
culations for the two neighboring values of x, 
in order to follow the effect of changing asym- 
metry on intensities in the final construction of 
the spectrum. 

This set is now sorted’? on J”, K’’_,, K’’; in 
the same order as the cards for the lower energy 
levels. The latter set is collated with the former. 
In this operation an energy-level card for a given 
value of J”, K’’_,, K"’; is placed in front of the 
group of cards in the transition set having the 
same values for J’, K’’_;, and K’’;. The com- 
bined set is now put through a Reproducer 
which gang punches information from the 
energy-level cards onto the transition cards, 
namely E’’(a’’, 6’, c’’) and other information 
to be described below. The transition set is now 
sorted on the quantum numbers of the upper 
state,? J’, K’_,, K's, and collated with the set 
for the upper energy levels, from which 
E’(a’, b’, c’) is gang punched. 

The next step is to gang punch*® onto the 
whole set the value of the band center, vo. The 
pack of cards is then put through the Multiplier, 
which carries out the addition and subtraction 
(1) to give the line position in one operation. 


Calculation of Line Intensity 


In order to calculate the line intensity the 
Boltzmann factor has first to be determined. The 
set for the lower energy levels is group multiplied 


8 At present there is no straightforward way to insert a 
quantity which is to be added to all cards directly into the 
machine to have the counterpart of group multiplying. If 
vo is not known beforehand, and only trial values are 
inserted in (1), we have found it convenient to avoid 
wasting columns for vo by assuming a center of 10,000 
cm™!, It is not necessary to punch this onto the cards, but 
merely to introduce unity into the fifth position of the 
summary counter by an impulse from the unit multiplier. 
The value of 10,000 rather than 0 was chosen to insure all 
calculated »’s are positive, whereby sorting is simplified. 


by 1/k7T. The cards are arranged in order of 
increasing E/kT on the sorter and collated with 
a set of cards representing the table of the 
exponential function, the latter being the pri- 
mary set. The value of exp —E/kT is gang 
punched onto the energy-level cards. 

The next step is to insert the appropriate 
weight factors gy,,. These depend on the 
symmetry classification of the energy levels. 
For instance, in the case of H2O the levels with 
7 odd have a weight of three, and the levels with 
r even have a weight of one. The cards can be 
put through the Sorter which separates the even 
and odd 7’s. The Boltzmann factors for the 
r-odd set are multiplied by three; in the r-even 
set the Boltzmann factor is multiplied by unity 
in order that the combined factor appears in the 
same field on all cards. The whole set is now 
sorted’ on J”, K’’_,, and K”’; for the collation 
with the transition set described above. The 
Boltzmann factor is copied from this set at the 
same time as the energy levels. 

For pure rotational spectra the quantity 
v(1—e-”/*7) appearing in (3) is important. It can 
be calculated on the transition set by the fol- 
lowing operations. First the exponential is 
evaluated by multiplving v of (1) by 1/kT, col- 
lating with the exponential table, gang punching 
e~"kT, Since the Multiplier can carry out the 
operation C—AB in one run, the required 
quantity can be evaluated with C=A=vy and 
B=e-"*T, all of these factors having been 
punched on the cards. 

‘ The line intensity can now be obtained on the 
transition set by putting through the Multiplier, 
which evaluates (3). This computation can be 
verified by sorting on the quantum numbers of 
the lower state and summing the intensities of 
all lines arising from each lower level. The sums 
should be equal to 


(2J+1) gs. exp —E”(a", b”, c’’)/kRT. (5) 


The product (5) is carried out on the lower- 
energy-level cards and listed. 

The final result is a set of cards bearing the 
quantum numbers, line positions, and line inten- 
sities for each transition, constituting all the 
lines of the spectrum up to® J<13. These cards 


®For molecules, such as H:O and HS, where one 
moment of inertia is small, the Boltzmann factor causes 
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Fic. 2. The a-band represents the 10,100A band of H2S (with a slight change in the intermediate 
inertia as described in footnote 11) as it would appear with a slit width of 1 cm™ at infinite resolution. 
(Band center 9911 cm~!.) The b-band respresents the appearance of a band under the same conditions 
arising from a transition between the same vibrational states if the moment were along the intermediate 
axis. (Band center arbitrarily taken as 10,000 cm~'.) Note that in the asymmetric rotor there is not a 
very marked difference between the two types of bands when observed with a finite slit width. 
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Fic. 3. Same bands as Fig. 2 with slit width 2 cm™! 


can be sorted and the information listed in a 
variety of ways. 


Ill. SPECTRUM ANALYSIS 
Completely Resolved Spectra 


By sorting on AJ, AK_;, AK,, and K_; or K,, 
the lines can be listed in branches, sub-branches, 
and wings (1I—Appendix), by means of which the 
structure of a band is analyzed. Inspection of the 
rate of decay of intensity in the wings will 
indicate whether any strong lines occur above 
J=12, and the intensities of such lines can often 
be found by extrapolation. 


all lines involving J>12 to be of negligible intensity. A 
satisfactory study of the bands of molecules with all 
moments large will require the extension of the tables of 
reduced energy levels and line strengths to higher J’s. We 
hope to do this by means of the Correspondence Principle. 


The cards can also be sorted with the lines 
arranged in order of increasing or decreasing 
wave number; or the strong lines, say with 
intensity greater than some predetermined value, 
can be listed separately. These printed lists are 
extremely convenient in analyzing well-resolved 
spectra. 


Unresolved Spectra 


In the analysis of incompletely resolved spec- 
tra, punched-card methods can be used to 
further advantage. A set of cards bearing con- 
secutive wave numbers in the region of interest 
is collated with the set of transition cards 
arranged in order of increasing wave number. 
The wave-number set has X punches in various 
columns so that when the combined set is put 
through the Tabulator, successive totals can be 
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Fic. 4. Same bands as Fig. 2 with slit width 10 cm~!. 


listed, giving the sum of the intensities over any 
interval of wave numbers corresponding to a 
slit width. Additional X punches are available 
to give the intensity at intermediate positions 
of the “slit.” A plot of the results gives the 
appearance of the band at the resolution of the 
slit width equal to the interval chosen. 


IV. TYPICAL RESULTS 


In order to test the punched-card procedure 
we have calculated the line positions and inten- 
sities J/<13 for the 10,100A band of H.S.'° The 


10P, C. Cross, Phys. Rev. 46, 536 (1934); 47, 7 (1935); 
J. Chem. Phys. 5, 370 (1937); P. C. Cross and B. L. Craw- 
ford, Jr., ibid. 5, 371 (1937); B. L. Crawford and P. C. 
Cross, ibid. 5, 621 (1937). 


constants chosen! were those given by Cross 
for the greatest and least moments, but the 
intermediate moments were adjusted to make 
x’ and x’ exactly 0.5 and 0.4, respectively. Cal- 
culations for x’=0.5 were done simultaneously. 
The plots in this paper refer to x’=0.4 only. 

A list of calculated positions and intensities, 
omitting all lines with intensities’? <1.0, was 


11 The values chosen for computation were 


a’ =10.393 a’=9,.910 

(b’’ =8.976) (b’=8.278) and (b’=8.550) 

c”’ =4.723 c’ =4.470 

x’ =0.5 x’ =0.4 and x«’=0.5. 
Cross’s value for b’’ was 9.040, giving «x’’ =0.523; and for 
b’ was 8.379, giving x’ =0.437. The value of T was 298°K. 
The band center was taken as 9911.03 cm™ for the cal- 
culations of the a-type band, and arbitrarily chosen as 
10,000 cm~ for the b-type. 

2 Actually, relative intensities I”, +;.’, ~/const. 

v(1—e-’/ET) were used. 
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compared with the observed lines.!° More than 
90 percent of the observed lines were identified 
and accounted for. The exceptions were lines 
found in a subsidiary list of those with intensities 
lying between 0.5 and 1.0. A few lines whose 
calculated intensity is ~3 were not given in the 
original list of observed lines. The agreement in 
wave number was satisfactory considering the 
rigid-rotor approximation of the calculations. 

The results of collating the wave-number set 
with progressive totals listed for successive unit 
wave-number intervals (slit width of 1 cm~) are 
plotted in Fig. 2. The curves for a resolution of 
2 cm™ are given in Fig. 3. The synthetic spectra 
obtained by this method correspond very closely 
in character to the experimental curves. For 
example, the 2.63u band of Fig. 1 (obtained with 
a slit width of 1.8 cm~') compares favorably 
with the b-band of Fig. 3, although the constants 
chosen in the latter are those for the 10,100A 
band. 

One immediate result of importance is that the 
observed band seems to have the appearance of 
the type where the electric moment is parallel 
to the intermediate axis, contrary to predictions 
based on symmetric-rotor formulas. 

By changing the control of the Summary 
Tabulator, sums were, obtained for an interval 
of ten wave numbers. The results are plotted in 
Fig. 4. The appearance of the bands is charac- 
teristic of the envelopes observed experimentally 
at this resolution. 


V. CONCLUSIONS 


The punched-card method: of obtaining line 
position and intensity, and of listing the sums 
of the intensities within any arbitrary slit width, 
gives curves which correspond very closely in 
character to band envelopes observed experi- 
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mentally at all degrees of resolution. This 
method remains accurate when the slit width is 
of the order of the average line separation, so 
that the detail of irregularities characteristic 
of the asymmetric rotor is preserved. 

The method also has the advantage over other 
procedures that the exact values of the position 
and intensity of the lines of the rigid rotor are 
used and are not averaged by replacing sums by 
integrals. This is particularly advantageous in 
parts of the spectrum involving levels of low J. 

With the procedure for drawing up spectra 
mechanized, dimensions of the molecules in the 
upper and lower states can be obtained, with an 
accuracy dependent on the detail of the band, 
by successive approximations, that is, by a 
process similar to that used in interpreting elec- 
tron-diffraction curves. 

Even when the resolution is poor, band centers 
can be located quite accurately. This is of value 
in fixing the vibrational levels of the molecules, 
precise values of which are of importance in cal- 
culating thermodynamic quantities from the 
partition function. 

The interpretation of the infra-red bands in 
terms of characteristic vibrations of the mole- 
cule is greatly dependent on the presence of 
certain branches or on the shape of the en- 
velopes. With a judicious choice of constants 
even the first run is sufficient to identify the 
direction of the electric moment with certainty 
and to locate the center of the band accurately. 
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N principle it is possible to determine the mo- may also happen, in consequence of the restric- 
lecular weight distribution of a high polymer, _ tion of the data to a certain concentration range, 
containing no excessively large molecules, from that components at the ends of the distribution 
interpretation of its light scattering in solutions curve will not be detected and, moreover, that 
at various concentrations. Actually there appear the conclusions regarding the remainder of the 
to be severe limitations to this method. Not only distribution will be distorted. 
is it evident that, considering experimental error, In particular it may be difficult or impossible 
the resolving power is low for the detection inany to perform a satisfactory extrapolation to infinite 
part of the molecular weight distribution curve dilution. In a recent letter to the Editor of this 
of a group of components of similar molecular Journal! and in a later article? Zimm, Doty, and 
weights (as contrasted to a single component), it Mark have used the expansion formula 
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Here H is a constant, 72 the turbidity due to the to deviate markedly from the parabola at concen- 
solute, c the concentration of the polymer in trations approaching one percent, and to be 
grams per milliliter of solution, f; the weight much less curved. It is therefore evident that 
fraction in the polymer of the species of molecular expansion (1) may be misleading even at reason- 
weight M;, and B and C are constants occurring ably low concentrations. 

in the expression for the osmotic pressure (see It is to be emphasized that expansion (1) may 
below). The weight average molecular weight is _ either be satisfactory at moderate concentrations, 
denoted by M, and ((AM/M)*)w and ((AM/M)*)» somewhat misleading as in Zimm, Doty, and 
are the second and third moments of the molecu-  Mark’s example, or wholly inadequate, depending 
lar weight distribution curve about the weight upon the nature of the molecular weight distri- 
average molecular weight. Expansion (t) may bution and osmotic behavior (e.g., the values of 
however fail even at the low concentrations used B, C, --+ in Eq. (1)) of the polymer. We now give 


4 H 1 . . . e 
by the ne authors. For oe rs “er “ some rather general considerations which possibly 
experimental points are plotted lor the light 244 something to our understanding of this 
scattering of -a mixture of two homogeneous situation 
polystyrene fractions and compared with the The equations which must be handled are 


parabola corresponding to the right-hand side considerably simplified and the considerations in 


of (1). If, however, the complete expression on jm ; 
the left-hand side of (1) is plotted for the mixture "°2*™ ste facilitated by the use of P. Debys's 


of the two fractions, the resulting curve is found 1B. H. Zimm and P. M. Doty, J. Chem. Phys. 12, 203 

















SS (1944). 
* Contribution No. 1024 from the Gates and Crellin ? P. M. Doty, B. H. Zimm, and H. Mark, J. Chem. Phys. 
Laboratories of Chemistry. 13, 159 (1945). 
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_ Fic. 1. Plot of L(y) for various two-component distribu- 

tions. Curve 1 represents a monodisperse system, curves 
2-5 (heavily drawn) tend toward the upper asymptote, and 
curves 6 and 7 (lightly drawn) belong to the lower 
asymptote. The limiting slopes are also indicated. The 
parameters involved are the following: 


M (()) Mi M: 
mu. \\it)/y Moh 
0 


(1) 1 1 0 1 
(2) 2 } 1.14 0.11 0.86 
(3) 2 ; 1.71 .29 50 
(4) 2 2 4.56 .44 .14 
(5) 2 16 32.5 .50 .016 
(6) 1.5 3 2.0 50 .333 
(7) 1.5 16 48.3 66 .007 


suggestion’ that Hc/r2 be considered as function 
of a quantity y(c), 


c 1 
T2 IiM; 
i 1+yM; 


where y is defined as follows. Let the osmotic 
pressure be given the representation 


=L(y), (2) 


fi 
r=RT|E “e+ o(0 ’ (3) 
in which g(c), whose first derivative vanishes 
at zero concentration, is independent of the 


’P. Debye, General Report No. 17, August 28, 1944. 
War Production Board, Office of Assistant Rubber Director 
of Research and Development of Synthetics. 
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degree of polymerization;‘ then y is given 
by y(c)=(d/de)y(c). (In (1) p=(2B/RT)c 
+(3C/RT)c?+--- was used.) The form of ¥(c) 
can be determined either by osmotic pressure 
studies or by measurements of the light scattering 
of monodisperse fractions. 

The following general remarks can be made 
about L(y): 

(a) The function (3) has always a positive 
slope 








aL fiM? 
— = 1? (4) 
ay i (it+yM,)? 
and is never curved upward, since 
eL Si:M:f;M; 
oe ae (M;—M;)?<0. 
oy? i (1+yM,)3(1+y¥M;)* 


(5) 


It is evident that in general 6°L/dy”, and hence 
the deviation of L(y) from a straight line, is 
greater, the greater the heterogeneity of the 
mixture. 

(b) For a polydisperse system L(y) is a curve 
whose initial slope is (1+((AM/M)*),)? and 
whose intercept is the weight average molecular 
weight. For the theoretical limit of large y, such 
that yM;>1, L(y) approaches an asymptote of 
unit slope, whose intercept is the reciprocal 
number average molecular weight 1/M,. Ac- 
cording to this and Eq. (5) the slope of L(y) is 
always greater than unity. 

(c) For a monodisperse system L(y) is simply 
a straight line of unit slope whose intercept is the 
reciprocal of the molecular weight. 

To show some of the wide possible variations 
in the form of L(y) we present theoretical plots 
for mixtures with various values of M/M, and 
involving solute molecules of only two different 
degrees of polymerization. The graphs are in 
terms of ML(y, M;), M.:/M and WM, since 
ML, M;)=L(My, M;/M). (Figure 1.) 

It should be noted that for any solution obeying 
Eq. (2) (and the curves given here for bimodal 
distributions illustrate this) any measurement of 
L provides an upper limit for M, and a lower 


‘For experimental evidence that this is a good approxi- 
mation cf. e.g., M. L. Huggins, J. Am. Chem. Soc. 64, 1712 
(1942); T. Alfrey, A. Bartovics, and H. Mark, J. Am. 
Chem. Soc. 65, 2319 (1943). 
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limit for M, while any two measurements of L at 
different y establish a lower limit for ((AM/M)*),. 
Whether these limits lie close to the actual values 
of the corresponding quantities depends on the 
curves and the values of My in question: for any 
range of My with ((AM/M)*)y small enough 
compared to M/M, satisfactory extrapolation to 
y=0 (as with the aid of Eq. (1)) is feasible for 
determining ((AM/M)*)4 and M, and with 
((AM/M)?)x,/(M/M,,) great enough the estimate 
of 1/M,, will be satisfactory, but generally the 
range of y in which significant measurements can 
be made may not be wide enough either to satisfy 
both these conditions or to make certain whether 


one of them is really satisfied. If extrapolation to 
yY=0 appears to be possible it may still be 
unsatisfactory because of the presence of a 
relatively small amount of material of relatively 
extreme high molecular weight such that the 
correct curve of L(y) bends significantly down- 
ward at values of ¥ lower than that of the lowest 
measurement. On the other hand, in the presence 
of material of relatively extreme low molecular 
weight, the curve of L(~) may have a maximum 
slope very little greater than unity and may seem 
to have reached an asymptote at the highest 
attained value of ¥ when actually its asymptote 
lies higher. 
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Theory of Multimolecular Adsorption from a 
Mixture of Gases 


TERRELL L. HILL 
Department of Chemistry, University of Rochester, Rochester, New York 
October 29, 1945 
HE purpose of the present note is to outline a possible 
way of extending the Brunauer-Emmett-Teller! 
(BET) theory of multimolecular adsorption to the case of 
a mixture of two gases. The problem will be discussed in 
more detail in a later paper, including applications and the 
treatment for any number of gases. 

We consider the adsorption of two gases A and B at 
relative pressures xa=pa/pa® and xp=pp/pp®, and at a 
temperature T below both critical temperatures. Let ga, gn, 
GAA, BB, QAB, YBA, 7ANa, and gpwa, be, respectively, the 
heats of adsorption per mole of A on the bare surface, B on 
the bare surface, A on a layer of pure A, B on a layer of 
pure B, A ona layer of pure B, B on a layer of pure A, A on 
a layer of composition Na (mole fraction of A), and Bona 
layer of composition Na. Let s; represent the surface area 
covered by i layers of adsorbed molecules. Then, at 
equilibrium, 


oro exp (—qa/RT), 
apsops = bpsi(1—Na,) exp (—gp/RT), 

oe basiNa; exp (—gawaci-»/RT), 

apsi-:pp = bpsi(1— Na,) exp (—gBvaG-»/RT). 

The above equations use average heats of adsorption and 
they assume that the energy of attachment of a molecule 
on top of the (¢— 1)th layer (to become a member of the ith 
layer) depends on the composition of the (¢—1)th layer 
only. If we consider similar equilibria for the evaporation 
and condensation of A from pure A and B from pure B, we 
have 


(1) 





MA _(Sitset:+++5n)Nait (set: ++ +5n)Nast:: 


+s,Na, 
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aaspa®=bas exp (—gaa/RT), 
apspp’=bgs exp (—gpp/RT) 
or 
b 
-A = pa° exp (gaa/RT), 
= (2) 


72. pp! exp (qpp/RT). 


We make the approximation, as is done in the BET 
theory, that second and higher layers have the same 
evaporation-condensation properties as the liquid state. 
Thus, in the present case, it seems appropriate to assume 
that gaw, and gpy, can be obtained from the vapor 
pressure curves for a binary liquid mixture of A and B at 
the same temperature. As a good approximation,” we can 
write for the equations of these vapor pressure curves 

*a= Na exp [(gaa—gas)(1—Na)*/RT], 
xp =(1—Na) exp [(gsp—gpa) Na?/RT]. 
One can of course use other equations in place of Eqs. (3), 
or, in fact, use the experimental vapor pressure curves if 
available. We can get corresponding equations involving 
gana and gpw, by considering equilibria of the type used 
to obtain Eq. (2): 
aaspa=basNa exp (—gan,/RT), 
apspp =bps(1— Na) exp (—gpy,/RT); 
or, from Eq. (2), 
*xa= Na exp [(gaa—gana)/RT], 


(3) 


4 

xp =(1—Na) exp [(@pp—gpwa)/RT]. @ 
Comparing Eqs. (3) and (4), we see that 

QAA—QANa= (gaa— as) (1—Na)?, (s) 


BB —@BNa= (¢BB—QBa) Na?. 


If Eqs. (2) and (5) are now substituted into Eqs. (1), one 
finds 


$¢=5;-1E;_1=SoEoEiE2:+-Ej-1, 1=1,2,--- (6) 
Eo=xaQa+xsQs, 
E;=xaQap? +xpQpae™, i= is ‘e — (7) 
[a(t)=Na?, B(t)=(1—Na,), 
Na, =% /Eo, ‘ 
Mam znOnaPOD (Ba, pene a 
Qa=exp [(ga—gaa)/RT], 
Qs =exp [(¢s—gss)/RT], 9) 


Qan=exp [(gan—gaa)/RT], 

Qpa=exp [(qBa—gsp)/RT]. 

If va and v,4 are the volume of A adsorbed and the 

volume of A necessary to give a complete unimolecular 
layer, respectively, then? 


(10) 





Sotsit s+ *+Sn 


Um* 


Fol (1+ 2it Ek 2+: ++ +EiEs: ++ En-1)Nait(Eit::++EiEs: + +En1)Nagt:+++Ei-++En-1Nay. 


(11) 





1+4£o(1+2£,+2,;E.+ ~ +FiF2: % *En-1) 


From Eqs. (7) and (8): 
lim Na;= Na=xaQan*/E, 


lim E;=E= xaQap?+xpQpa“. 


to 


(12) 
(13) 





It will be noted that vA/v,4 and vB/v,® become infinite 
(condensation) when E=1 and m=, Putting E=1 in 
Eqs. (12) and (13) gives Eqs. (3), as expected. 

It appears that Eq. (11) may be approximated to within 
several percent in most cases of interest by 
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Eo{ Na, (E*!—E*—E+1) 
~~. + Nal (n—1)E**!—nE*+E}} (14) 
Tn (1—E£)[1—E—E,(E*—1)] 
and, for n= «, by 
va Eo Nas(1—E)+NaE}, - 





Ym’  (1—E)(1— E+ Eo) 


The quantity vB/v,,® may be obtained from Eqs. (10), 
(11), (14), and (15) by replacing Na by 1— Na. 

If xp=0, these equations reduce to the usual BET 
equations with Na=1, E=xqa=x and Eo=xaQa=xc. 

In comparing theory with experiment, Qa, Qs, ¥m4, and 
0m® may be found from the pure A and pure B isotherms, 
while the vapor pressure curves can be used to give Qap 
and Qpa. 


nn P. H. Emmett, and E. Teller, J. Am. Chem. Soc. 60, 309 
2 See, for example, R. H. Fowler and E. A. Guggenheim, Statistical 
Thermodynamics (Cambridge University Press, London, 1939), p. 357. 
3 It is an approximation to use mole fractions instead of area fractions 
in Eqs. (10), (11), (14), and (15) unless the molar volumes of A and B 
are equal. This will be discussed in a later paper. 





A General Theory of the Reaction Loci in 
Emulsion Polymerization. II* 


WILLIAM D. HARKINS 
Department of Chemistry, University of Chicago, Chicago, Illinois 
October 18, 1945 


HE earlier theory! has been extended to explain how 

two minor loci can act, either after soap micelles 

disappear in polymerization with soap or when no initial 
soap is used: 


A. Loci for the initiation of polymer particle nuclei. 

1. Soap micelles: When an initial soap concentration 
of about 3 percent is used (with 100 g monomer to 
180 g water), nearly all of the polymer particle 
nuclei are initiated in the micelles. With more soap 
the polymer particles at any given yield are more 
numerous and smaller, and with less soap less 
numerous and larger. For example, an isoprene- 
styrene polymer at 40-percent yield gave diameters 
of 375A and 485A in 4-percent and 1-percent 
potassium myristate solutions. 

2. Water: The aqueous phase, excluding soap micelles 
and polymer-monomer particles, is a secondary 
locus whose relative prominence increases as soap 
disappears from its micelles due to adsorption of 
soap molecules around the monomer-polymer 
(latex) particles. Finally (e.g., at a 60-percent 
yield) by process 4, below, too few molecules of 
monomer remain to make this locus effective. 

3. Monomer droplets: The new evidence indicates that 
monomer droplets are relatively unimportant in 
initiating new polymer particles, since thus far no 
evidence has been obtained to show that a monomer 
droplet can do more than change into a polymer 
particle of much smaller size. 
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B. Loci for the formation of polymer. 


4. Polymer-monomer particle: Almost all of the 
polymer is produced in this locus. (See Com- 
munication I.) 

5. Monomer-emulsion droplets: The principal function 
of the monomer droplet is to act as a reservoir from 
which monomer diffuses through a seemingly thick 
diffusion layer into the other loci. In-static experi- 
ments every emulsion monomer droplet observed 
changed into a very much smaller polymer droplet 
when suspended in a soap solution at 40° to 50°C 
with all of the necessary reagents present. 


Water: Evidence that water (with 0.2 percent of po- 
tassium persulfate as catalyst) can act as a locus for the 
initiation of polymer particle nuclei is given by experiments 
(1) by Stearns and (2) by Corrin. (1) A flask partly filled 
with this aqueous phase has two side tubes, turned down- 
ward, sealed into the neck. In the bottom of one of these is 
isoprene, or styrene, or other monomers, and in the bottom 
of the other a mercaptan. The system is evacuated through 
the neck of the flask, which is then sealed off, after which 
the whole system is suspended inside the water of a 
thermostat held at 50°C for 18 hours. In two experiments 
of this type with styrene, the diameters of the polymer 
particles obtained by Oppenheimer were 1500 and 2000A, 
respectively, while the numbers were 37 X 10" and 5.6 X 10" 
per cu. cm: with isoprene the values were 2200A and 
2.910". (2) The polymerization of styrene in water and 
salt solutions was investigated in the absence of the water- 
monomer interface. Tubes with sealed-in inner cups were 
employed to insure a short vapor path between the bulk 
styrene and water. The tubes were maintained in a 40°C 
bath for 42 hours. The solutions contained 0.1 percent 
potassium persulfate, and any salt added was potassium 
chloride. (Table I.) 

Since at 50°C a molecule of mercaptan in the water is 
readily converted by the catalyst into a free radical (in (2) 
no mercaptan is present, so styrene itself forms the initial 
free radical), which can unite with a monomer molecule to 
give a larger free radical and again with more monomer 
molecules, until the growing chain is terminated, the 
results cited above accord with expectations. They indicate 
that the aqueous phase, after the disappearance of soap, 
can act as a locus for initiation of polymer particles. 

Locus of polymerization in monomer droplets (Stearns). A 
drop of styrene or isoprene with 0.2 percent lauryl mercaptan 
was suspended in an aqueous solution of 3 percent potas- 
sium laurate, 0.3 percent potassium ferricyanide, 0.1 
percent potassium persulfate, and 1 percent sodium car- 











TABLE I. 
Mean 
» X10? particle Number of 
% (u=ionic Conversion diameter particles 
K2S20s8 m KCl strength) (mg/cc) inA per cc 
0.1 0.0 1.11 2.89 1710 1.1 X10" 
0.1 0.0025 1.36 2.99 1310 2.5 X10" 
0.1 0.005 1.61 3.08 1370 2.3 X10" 
0.1 0.01 2.11 2.82 2870 2.3 X10" 
0.1 0.05 6.11 2.47 * * 
0.1 0.10 11.1 1.74 * * 
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bonate at 40° or 50°C. The principal phenomenon observed 
was that 0.999 of the volume disappeared by diffusion of 
the monomer into the soap solution. However, in every one 
of numerous experiments the monomer droplet gradually 
changed into a monomer-polymer droplet. 

After the drop was reduced to about a tenth its original 
diameter, diffusion into the soap solution ceased, since the 
droplet remaining consisted largely of polymer. In some 
cases the density was higher than that of water (plus soap), 
shown by its falling downward through the water. In gen- 
eral, the drop became deformed, becoming non-spherical 
as its content of polymer increased. 

* The work reported in this paper was done in connection with the 
Government Research Program on synthetic rubber under contract 


with the Office of Rubber Reserve, Reconstruction Finance Corporation. 
1W. D. Harkins, J. Chem. Phys. 13, 381 (1945). 





The Density and Thermal Expansion 
of Liquid Phosgene 


C. N. Davies 


Industrial Health Research Board, Medical Research Council, 
London, England 


December 10, 1945 
1. INTRODUCTION 


OME unpublished values for the density of liquid 

phosgene suggested that the figures available in 
chemical literature might be high. The measurements 
described in the present paper confirmed this view; they 
were carried out on two samples, from the same source, 
which had been shown by chemical analysis to contain 
only traces of impurity too small to affect the density 
within the accuracy of the determinations. 


2, EXPERIMENTAL METHOD 


A bulb was blown on thick-walled Pyrex capillary tubing 
and the volume per cm of stem determined over the work- 
ing range by using a cathetometer to measure distances of 
the meniscus of mercury charges of known weight from a 
mark on the tube. The volume up to this mark was 1.7505 











TABLE I. 
First series * Second series 
Temperature ensity Temperature Density 
—9.35°C 1.444; g/cm? 20.05°C 1.3730 g/cm* 
—8.60 1,442) 20.10 1.3714 
—7.75 1.4397 27.00 1.3547 
—5.70 1.4348 37.00 1.3310 
—3.25 1.429 45.45 1.309% 
—2.70 1.4277 48.50 1.3019 
—2.00 1.4260 19.20 1.3745 
—0.65 1.4224 14.50 1.3853 
—0.05 1.4214 5.00 1.4076 
15.03 1.3839 0.00 1.4187 
16.73 1.3803 —5.00 1.4285 
20.59 1.3709 —10.00 1.4403 
28.01 1.3533 —15.00 1.4513 
34.66 1.3362 —20.00 1.4631 
40.13 1.3229 20.00 1.3714 
45.27 1.3093 
50.75 1.2951 
55.50 1.2831 
60.20 1.2706 
17.75 1.3783 
12.00 1.3914 
7.55 1.4022 
—0.20 1.4207 
—6.10 1.4348 
—10.40 1.4450 
—15.50 1.4576 
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+0.0004 cm* and the tube diameter was about 4 mm. 
Phosgene was distilled in and the tube sealed. The posi- 
tion of the phosgene meniscus was then determined at 
temperatures from —20°C to +60°C. The mass of phos- 
gene was corrected for the mass in the vapor phase using 
the vapor pressures of Paternd and Mazzucchelli.! Correc- 
tions were also applied for the meniscus, and weighings 
were reduced to vacuum. Two series of experiments, in the 
order given below, were carried out independently on 
different samples. 

Possible errors caused by the expansion of the bulb 
under pressure and temperature were shown to be negligi- 
- 3. EXPERIMENTAL RESULTS 

The results given in Table I, taking both series together, 
have been fitted by the method of least squares with the 
following regression curve: 


1.42014 —0.0023120¢ —0.00000287 22? 
= 1.42014 (1 —0.001628¢ —0.000002022/*), 


where p is the density in gm/cm at £°C. 

By differentiating this equation the thermal expansion 
per degree centigrade was found to have the following 
values: 


p 


Temperature 1/0(do/dt) 
50°C 0.002003 

30 0.001842 

10 0.001696 

—10 0.001563 


Above 0°C the agreement between the two series and of 
each with the equation is better than 0.1 percent. Below 
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zero the scatter increases but is always less than 0.2 per- 
cent. This may be attributed to more accurate temperature 
control during the determinations at the higher tempera- 
tures. 


4. COMPARISON WITH PREVIOUS DETERMINATIONS 


Paternd and Mazzucchelli! give an equation based on 
experiments over the range —15.4°C to +59.9°C which is 
about 0.6 percent above the present data at 20 C and 
gives a slightly lower coefficient of expansion. 

Atkinson, Heycock, and Pope? found still higher den- 
sities and exceed the latest figures by about 1.2 percent. 

Germann’ only did determinations at 0°C and 25°C and 
agrees closely with Paternd and Mazzucchelli. 

Figures given by Emmerling and Lengyel‘ and Beck- 
mann’ are fairly close to Atkinson, Heycock, and Pope. 

All these determinations may be compared with the 
present findings (see Fig. 1). It is significant that the 
older values are high which suggests that the presence of 
carbon tetrachloride might explain the discrepancies. 
About 6 percent would raise the density of phosgene by 1 
percent and the deviation would tend to increase with rise 
of temperature as is observed in the case of the Italian 
data. It is probable that dissolved chlorine would also 
raise the density. The results of Atkinson, Heycock, and 
Pope were obtained by chemical estimation of the phos- 
gene instead of direct weighing. It is not possible to draw 
any conclusions from the methods of manufacture of the 
phosgene since various processes have been used and no 
correlation with the observed densities is apparent. 


5. SUMMARY 
A determination of the density and thermal expansion 
of liquid phosgene over the range —20°C to +60°C is 
described. The density in this range is given by the follow- 
ing equation: 
p = 1.42014 —0.0023120¢ —0.00000287 22%. 
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Statistical Mechanics of Linear 
Association Equilibria 


PAuL J. FLtory 


Research Laboratory, The Goodyear Tire and Rubber Company, 
Akron, Ohio 


November 21, 1945 


N a recently published paper Tobolsky and Blatz! 
have discussed the “possible application’’ of the 
lattice model for mixtures of linear polymers with small 
molecules to reversible associations of liquids such as 
the alcohols and their solutions in non-associated solvents. 
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It appears that they were unaware that this particular 
problem can be treated by previously derived relationships 
pertaining to heterogeneous polymer-solvent mixtures.* 
Their paper contains the unfounded and misleading state- 
ment, “ ... whereas Flory? was considering stable 
polymer molecules which did not change their size distri- 
bution when diluted, we are here concerned with the 
problem of reversible association. . . . "’ It is the writer’s 
purpose to emphasize here the generality of the previously 
derived thermodynamic functions as applied to equilibria 
involving transformations between various polymeric 
species. 

In applying the methods of thermodynamics, or sta- 
tistical mechanics, to chemical equilibrium problems the 
character of the bonds, or intermolecular forces, subject to 
equilibration need not be specified. The treatment, previ- 
ously presented, of equilibria involving polymer molecules 
is completely general, and although secondary valence 
polymers were not mentioned, it should have been clear 
that it is also applicable to these as well, provided that the 
“polymers” consist predominantly of linear sequences of 
structural units. This treatment led to the conclusion that 
the equilibrium constant for processes involving polymeric 
substances assumes the same form as that for the analogous 
reaction between monomeric species; the relationship of the 
equilibrium constant to the standard state free energy 
change is the same in the two cases. Thus, no modifications 
in the formulation of equilibrium constants are occasioned 
by the configurational randomness of the polymeric reacting 
species. 

For the purpose of applying this principle to associative 
equilibria specifically, consider a solution containing No 
associating units per unit volume, each of which is capable 
of entering into two “bonds’”’ to form linear polymers. If 
each of these bonds forms independently of the presence or 
absence of other bonds attached to the adjoining units, 
then an equilibrium constant governing the associative 
process can be written 


K=p/(1—p)*No, (1) 


where p (equivalent to the “y” of Tobolsky and Blatz) is 
the “degree of association,” i.e., Nop represents the concen- 
tration of association bonds at equilibrium. Also from the 
above assumptions and the nature of the dependence of 
partial molal free energies on concentrations, it follows that 
the molecular size distribution is given by 


Nz=No(1—pyYp?", (2)4 
where N, represents the molar concentration of x-mers; the 
number average degree of polymerization, according to (2), 
is given by 

Z,=1/(1—), (3) 
the partial molal free energy of the solvent in this system is 
obtained directly by substituting (3) into Eq. (17) of 
reference 2, giving 

AF, =RT[In (1—v2)+02p +02" ], (4) 
where v2 is the volume fraction of solute including all 
species, and y is a semi-empirical parameter®® in which is 
included the heat of dilution term. Equation (4) is identical 
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with Tobolsky and Blatz’s Eq. (5), which is merely a 
special case of the more general equation, namely that in 
which the molecular size distribution of the solute is given 
by (2). Similarly, if desired, partial molal free energies for 
any of the various solute species can be derived by ap- 
propriately substituting Eqs. (2) and (3) into Eqs. (18) and 
(19) of reference 2. The resulting expressions differ from 
Tobolsky and Blatz’s (6) and (7), which latter contain the 
intractable quantity y’ and refer to the ambiguous standard 
state consisting of ‘‘pure mixed”’ solute molecules. 
Whereas Tobolsky and Blatz prefer to include the factors 
controlling ~, the extent of association (or reaction), in the 
partition function, the writer has preferred to relate p 
separately to an equilibrium constant. The procedures are 
fundamentally equivalent, the writer’s being generally more 
easily applied. 
The factor b"0+2n, introduced into their intuitively 
derived partition function to account for mixing of mole- 
cules with holes in the lattice appears to be both unneces- 
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sary and incorrect. Tobolsky and Blatz! include the 
entropy of bond formation (& In a[2(x—1)Nz]), which had 
been omitted by Tobolsky’ previously. Hence, their 
definition of P as the heat of association is correct, but in 
the previous paper’ P must be considered to represent a 
standard state free energy change, as previously pointed 
out.? 

In conclusion it should be emphasized that the previously 
derived thermodynamic relationships are general, and a 
new derivation of the partition function for each applica- 
tion to a specific problem would seem to be unnecessary. 


1A. V. Tobolsky and P. J. Blatz, J. Chem. Phys. 13, 379 (1945). 

2P, J. Flory, J. Chem. Phys. 12, 425 (1944). 

3P, J. Flory, J. Am. Chem. Soc. 58, 1877 (1936) ; ibid. 64, 2205 (1942). 

4K. H. Meyerand A. van der Wyk, Helv. Chim. Acta 20, 1321 (1937), 
derived Eq. (2) for association processes in which the free energy or 
for each successive monomer added to the micell is assumed to be the 
same. 

5 M. L. Huggins, J. Phys. Chem. 46, 15 (1942); Ann. N. Y. Acad. Sci. 
43, 1 (1942); sbid. 44, 431 (1943); J. Am. Chem. Soc. 64, 1712 (1942). 

6 Pp, J. Flory, J. Chem. Phys. 13, 453 (1945). 

7A. V. Tobolsky, J. Chem. Phys. 12, 402 (1944). 
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